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ABSTRACT: We review the current status of flavor-changing neutral currents in the charm sector. We 
focus on the standard-model predictions and identify the main sources of theoretical uncertainties in both 
D — D° mixing and rare charm decays. The potential of these observables for constraining short-distance 
physics in the standard model and its extensions is compromised by the presence of large nonperturbative 
effects. We examine the possible discovery windows in which short-distance physics can be tested and study 
the effects of various extensions of the standard model. The current experimental situation and future prospects 
are reviewed. 
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1 INTRODUCTION 

The remarkable success of the standard model in describing all experimental information to 
date suggests that the quest for deviations from it should be directed either at higher energy 
scales or at small effects in low-energy observables. To the last group belong measurements 
(with precision surpassing 1%) of electroweak observables at CERN's Large Electron Positron 
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Collider (LEP) and the SLAC Linear Detector (SLD) as well as the Tevatron experiments (1). 
Tests of the standard model through quantum corrections are a powerful tool for probing the 
high energy scales possibly related to electroweak symmetry breaking and the flavor problem. 
The absence of flavor-changing neutral currents (FCNC) at tree level in the standard model 
implies that processes involving these currents are a primary test of the quantum structure of 
the theory. The study of FCNC has been focused on processes involving K and B mesons, such 
as K°-K° and B^-B^ mixing, and on rare decays involving transitions such as s — > M + l—, 
s — ► dvv, b — ► S7, and b — > si + l~ . 

The analogous FCNC processes in the charm sector have received considerably less scrutiny. 
This is perhaps because, on general grounds, the standard-model expectations are very small 
for both D°-D° mixing and FCNC decays. For instance, no large nondecoupling effects arise 
from a heavy fermion in the leading one-loop contributions. This is in sharp contrast with 
K and B FCNC processes, which are affected by the presence of the top quark in loops. In 
the standard model, D-meson FCNC transitions involve the rather light down-quark sector, 
which implies an efficient Glashow-Iliopoulos-Maiani (GIM) cancellation. In many cases, 
extensions of the standard model may upset this suppression and contributions may be orders 
of magnitude larger than the standard-model short-distance contributions. 

As a first step, and in order to ponder the existence of a window for the observation of 
new physics in a given observable in charm processes, we must compute the standard-model 
contributions to such quantities. Attention to standard-model backgrounds is particularly 
important in this case owing to the presence of potentially large long-distance contributions 
that are nonperturbative in essence and therefore noncalculable by analytical methods. In 
general, the flavor structure of charm FCNC favors the propagation of light-quark intermediate 
states. Thus, if it turns out that the charm-quark mass is not heavy enough compared to a 
typical scale of hadronic effects, long-distance effects are likely to dominate. They will obscure 
the more interesting short-distance contributions that are the true test of the standard model. 
Large long-distance effects are expected in D°-D° mixing and FCNC charm decays. In the 
case of mixing, although the long-distance effects dominate over the standard-model short- 
distance contributions, there could still be a significant window between these and the current 
experimental limits. The predictions of numerous extensions of the standard model lie in this 
window. We examine this possibility in Section 2, where we also show that even the current 
experimental bounds on D° — D° mixing constrain several new-physics scenarios. 

Charm radiative decays are completely dominated by nonperturbative physics and do not 
constitute a suitable test of the short-distance structure of the standard model or its exten- 
sions. However, semileptonic modes such as c — > u£ + i~ may be used to constrain various 
standard-model extensions, since their kinematics might allow measurements away from the 
resonance-dominated region, where the bulk of the long-distance contributions lie. We demon- 
strate this in detail in the case of supersymmetric (SUSY) theories with or without imparity 
conservation. Purely leptonic flavor- violating modes such as D° — ► /i + e~ are free of standard- 
model backgrounds. We review the status and prospects of rare charm decays in Section 3. 



The time evolution of the D° — D° system is described by the Schrodinger equation as 
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where P denotes principal value, we see that the absorptive part, T, comes from summing 
over real intermediate states 
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On the other hand, the contributions to Mij will include not only the second term in Equa- 
tion 2 but also a dispersive piece from off-shell intermediate states. Because M and T are 
Hermitian matrices, M\ 2 = M| x and Ti2 = T^. Invariance under CPT requires M\\ = M22 
and Tn = T22- The Hamiltonian eigenstates in terms of the weak eigenstates are 
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Defining the "right-sign" amplitudes Aj= A(D° — > /) and Af = A(D° — > /), we can also 
define the "wrong-sign" amplitudes Aj = A(D° — > /) and Af = — ► /), where / and / 
are CP-conjugate final states. Then the time evolution of states that start as weak eigenstates 
at t = results in the time-dependent wrong-sign rate given by 
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The eigenvalues of the Hamiltonian are 71 2 = 
define 
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Section 2.4 discusses various ways of accessing the mixing parameters x and y. 
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1 Standard- Model Predictions for Mixing Parameter 

In this section, we review our understanding of the standard-model predictions for D° — 
D° mixing. As always in charm FCNC processes, the main issue is to estimate as accurately 
as possible the size of the long-distance contributions. Only then we can evaluate the potential 
of D° — D° mixing to test the standard model and constrain its extensions. 
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2.1.1 SHORT-DISTANCE CONTRIBUTION TO Am D 

In the standard model, the short-distance AC = 2 transition occurs via box diagrams. The 
effective interactions at the m c scale are described by the Hamiltonian (2): 

( 2 2\2 

cfr - ~4^2 |KcsKcd| ^2 (<2 + 20), (10) 

with O = ^7^(1 — 75)0 M7 M (1 — 75)0 and C = u(l + 75)0 u(l + 75)0. The additional operator 
O' arises from the nonnegligible external momentum. The matrix elements of the operators 
can be parameterized by 

{D»\Om = \mlf D B D ■ (D°\O'\D ) = -l^ 2 rnlf D B' D . (11) 

In the vacuum-insertion approximation, one has Bd = B' D = 1 and the box diagrams' contri- 
bution to the mass difference is 

Am b D ox ~ 1.4 x 10~ 18 GeV ( -J^—) ( f ° V. (12) 
D VO.lGeVy V0.2GeVy V ; 

The 6-quark contribution is much smaller owing to additional Cabibbo-Kobayashi-Maskawa 
(CKM) suppression. Thus, for typical values of fjj and m s , the box diagrams contribute with 
approximately a^box — 

(few x 10~ 6 ) - (10 -5 ). In contrast to K and B mixing, the internal 
quarks in the box diagrams here are down-type quarks. The 6-quark contribution, which 
would give in principle the largest GIM violation, is suppressed by small CKM mixing factors. 
The leading contribution, as shown in Equation 10, is given by the strange quark and therefore 
results in a very effective GIM suppression. 

The box-diagram contributions to Ar are further suppressed by m 2 s . This can be seen as 
the two powers of the helicity suppression factor from the matrix element of a 0~ meson to a 
pair of quarks. 

2.1.2 LONG-DISTANCE CONTRIBUTIONS 

The long-distance contributions to D°—D° mixing are inherently nonperturbative and cannot 
be calculated from first principles. It is however extremely important to estimate their size 
in order to understand the origin of a possible experimental observation. 

A first step is to recall that D° — D° mixing is an SU(3)-breaking effect (3, 7): 

x, y ~ sin6»2 x [SU(3) breaking]. (13) 

Thus, the task is to estimate the amount of SU(3) breaking. An important observation (4, 11) 
is that the SU(3) breaking is a second-order effect in the quark masses. This circumstance 
would lead to the naive estimate 

x, y ~ sin# c 2 x (^-) 2 <> O(10~ 3 ), (14) 

V A hadr. / 

with Ahadr. ~ 0(1) GeV a typical hadronic scale. There are two basic methods to estimate 
the contributions to mixing beyond box diagrams: the heavy-quark effective theory (HQET) 
approach pioneered by Georgi (4), which is essentially rooted in the operator product expan- 
sion (OPE); and an exclusive approach computing the contributions of hadrons in complete 
SU(3) multiplets to the dispersive and absorptive parts. 
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D° — D° Mixing in Heavy-Quark Effective Theory: The applicability of the HQET 
ideas to D° — D° mixing rests on the assumption that the charm-quark mass is much larger 
than the typical scale of the strong interactions. Georgi pointed out (4) that in this case there 
are no nonleptonic transitions to leading order in the effective theory, since they would require 
a large momentum transferred from the heavy quark to the light degrees of freedom. This 
means that, in the effective low-energy theory, mixing is a consequence of matching the full 
AC = 2 theory at the scale m c with the HQET and then running down to hadronic scales 
(<S m c ). No new operators enter at low energy. The "long-distance" effects are provided by 
the tree-level matching of operators that are nonleading in the l/m c expansion but that are 
less GIM-suppressed than the box diagrams. 

First let us consider the four-quark operators generated from the box diagrams by inte- 
grating out the Ws. The contribution of these operators to the mass difference behaves like 
(4) 

. (4) 1 mi . , 

Am ° - 16^4' (15) 

where the first factor comes from the loop and m^ is neglected. This presents the same amount 
of GIM suppression as the box diagrams and therefore causes no enhancement. 

There will also be higher-dimension operators that, although suppressed by additional pow- 
ers of l/m c , can give important contributions, can give important contributions. For instance, 
six-quark operators arise by "cutting" one of the light-quark lines in the loop and then shrink- 
ing the connecting line leftover. As a consequence, two powers of m s are lost and the contri- 
bution from six-quark operators goes like 

Amg>* — 4 k/ 2 ), (16) 
m c m~ V / 

where the last factor comes from taking the hadronic matrix elements and / is the pseudo- 
Goldstone-boson decay constant. This constitutes an enhancement of (47r/) 2 /m c m s with 
respect to the four-quark operator contribution. 

Finally, eight-quark operators are obtained by cutting the remaining light-quark line and 
bridging the two four-quark pieces with a gluon. The resulting contribution goes like 

Am? - (17) 
47r 

This is the least GIM-suppressed contribution, although it is suppressed by and, most 

important, by the factor a s /4-7r. Therefore, no enhancement is expected from these operators. 

Ohl et al. performed a very detailed calculation in this approach, including QCD corrections 
to one loop (5). Their results can be summarized as 



Am? ~ (0.8 - 1.5) x 10~ 16 GeV 



m. 
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Am? ~ (0.6 - 1.6) x 10~ 16 GeV ( ™ s V (18) 

\ O.lGeV / 

A?n? ~ (0.05 - 0.25) x 10~ 16 GeV ( - j^v f ■ 

The lower numbers correspond to adding the various coefficients in quadrature, since naive 
dimensional analysis does not predict the relative signs. The upper numbers come from adding 
all contributions coherently. This leads to 

x ~ (0.6 - 2) x 10" 4 (19) 
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for m s = 0.1 GeV. Thus, the "long-distance" enhancement observed in this approach is of 
about one order of magnitude. 

More recently, Bigi & Uraltsev made similar arguments also in the context of the OPE (6). 
However, they found a larger eight-quark operator contribution than Equation 18 indicates. 
This results in x ~ 10~ 3 . Bigi & Uraltsev also undertook (6) a study of Ar in the OPE 
approach. They found that the most important contribution to y comes from the eight-quark 
operators, where the gluon connecting the two sets of four-quark operators is split to generate 
an imaginary part in the diagram. The "suppression" to be paid (relative to Amp) is given 
by Pos/Att ~ 0(1), where (5 = (11 — 2ny/3) is the QCD beta function. This results in 
y ~ x ~ 10~ 3 . 



Exclusive Approach: An estimate of the long-distance contributions can be obtained by 
assuming they come from the propagation of hadronic states to which both D° and D° can 
decay. There will be one-, two-, three- . . . particle intermediate states. Each of these groups 
can be further separated into sets whose contributions vanish separately in the SU(3) limit. 
For instance, one of these sets is formed by the two-charged-pseudoscalar intermediate states 1 
7r + 7r~, K + K~ , K~n + , and K + ir~ . Thus, computing their contribution to the mass difference, 
as shown schematically in Figure 1, gives a concrete realization of the estimate in Equation 13 
for an SU(3) set for which data are available. This was first done (7) in the massless limit. 
Although these "self-energy" diagrams will depend on the interaction chosen for the vertices, 
they have a universal imaginary part that typically comes from a logarithm. 

The contribution of Figure 1 to the mass difference obeys a dispersion relation of the form (7) 



1 r°° im[E(s)] ds 
7T J s (s -p 2 + ie) 



n P 2 ) = - 1 ;: l t."° . (2°) 



where sq = (mi + TO2) 2 , and mi and m2 are the masses in the loop. The mass difference 
is then given by Amj) = —Re\^(m 2 D )\/2m£). Taking into account a subtraction forcing the 
condition S(0) = 0, and keeping the masses, the implementation of a cutoff A in the dispersive 
integral (Equation 20) gives (8) 

m D \ B(tt + tt-) t . a . B(K + K~) r . 
x ~ — — < l{mir, m^, A) H I(mK,m K ,A) 

47T I Pt,-,,- PKK 



rJB(K-ir+) B(K+tt-) r . .A , . 

-2 cos 5^- — — i mjf , A) \ , (21) 

PA'tt J 

where Pij is the magnitude of the three-momentum in the actual decay, 5 is the strong phase 
between the Cabibbo-allowed and the doubly Cabibbo-suppressed decays, B is the branching- 
fraction, and the integrals are 



A-^ds 

I(mi,m 2 ,A) = - I 1 (22) 

J Sq S TY^D 

By taking the massless limit in Equation 21, one recovers the result in Reference (7) with 
the identification /x 2 = 2mo (A-m^). Although the result depends strongly on the cutoff 
A, this can be interpreted as the value of s for which the internal momentum reaches its 
maximum. Not surprisingly, the value of A that gives an internal momentum of ~ 1 GeV also 
yields fx ~ 1 GeV. This is A ~ (2-2.2) GeV, not far above mp. Using this cutoff results in a 

1 Strictly speaking, this is a U-spin set. 
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contribution to the mass difference of 

x ~ 6.5 x 1(T 4 x [1 - cos 5 b], (23) 

where we define b by 



V cd V us 



v*v. 



ad 



B(K-ir+ 

as a measure of the amount of SU(3) breaking in the absolute value of the amplitudes, in 
addition to the SU(3)-breaking relative phase 5. Measurements in doubly Cabibbo-suppressed 
(DCS) decays result in b ~ 1.22 (9). 

It is clear that the cancellation in this SU(3) set is quite effective. Although the dependence 
on the cutoff A appears unsatisfactory, it has a clear physical interpretation. In order to 
derive Equation 21, we assumed that the couplings at the vertices in Figure 1 are pointlike 
and therefore constant in s. However, the hadrons participating in these interactions are 
composite, with the compositeness scale Ahadr. ~ Oil) GeV, a typical hadronic scale. Thus, 
at or near Ahadr.) the couplings should develop a form- factor suppression determined by this 
energy scale that in turn would render the integrals in Equation 22 finite without ad hoc 
cutoff. Then we should replace 7(mi,m2,A) with 



oo w 1 — ^ ds 

I( mi ,m 2 ) = -I g 2 (s)' =g— , (25) 

s — m , 



ls o llb D 

with the coupling of the interaction vertex normalized so that g{s) ~ 1 for s < m 2 D and it 
falls off with Ah a dr. f° r s > Of course, the dependence on the cutoff is now hidden in 

g(s), but its physical meaning is more transparent. 

In principle, pseudoscalar- vector (PV), vector- vector (VV), and all other possible interme- 
diate states can be treated similarly. However, in many cases, the data are incomplete. In 
some cases, phase space suppresses some of the modes, but unlike in Ar (11) (see discussion 
below), there is no SU(3)-breaking effect induced by running out of phase space. Modes that 
are not accessible to the actual D decay (and thus to the absorptive part) still contribute 
to the dispersive part. We cannot express these in terms of branching ratios because the 
decays are not physical. But in principle the amplitudes for these "off-shell" contributions 
are well-defined. We then conclude that estimates provided by lighter, far-from-threshold 
sets, such as the one in Equation 21, are likely to be the most reliable ones in this approach. 
The underlying assumption is that the contributions of other sets, their signs being not fixed, 
should not conspire to add coherently and change significantly the order of magnitude of the 
effect. This reasoning results in 

x & 10" 3 , (26) 

with large uncertainties. In sum, the exclusive approach seems to be consistent with our naive 
estimate of Equation 14. Additional sources of enhancement in D° — D° mixing are discussed 
in Reference (10). 

The exclusive approach can also be applied to ATp. For instance, for the PP (two pseu- 
doscalars) set considered above, the contribution to y has an expression similar to Equation 21, 
so that the cancellation is still quite efficient (11). However, y receives contributions only from 
real intermediate states, so phase space could be a considerable source of SU(3) breaking, in a 
way that did not arise for the x case (11). In the most extreme cases, some states belonging to 
a given set would simply be absent because Ds are not allowed to decay into them. This could 
upset the SU(3) cancellations, perhaps to a large extent. This source of SU(3) violation is 
estimated by considering the effect of phase space in the contribution of a given intermediate 
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state F and within a given SU(3) representation R. Thus, one can write 

= Y: n (D°\H w \n)p n (n\H w \D } = j: n (D \H w \n)p n (n\H w \D } 
VF,R En(D°\H w \n)p n (n\H w \D^ En^(D° - n) ' [ ' 

where the sums are performed over intermediate states n belonging to the class F (PP, PV, 
etc.) and the representation R. Here, p n is the phase space of state n and depends only on 
the masses of the particles in this intermediate state. The expression in Equation 27 would 
be y if the states n of (F, R) were the only states available for the D° to decay into. Then, 
in order to obtain the corresponding contribution to y, yF,R must be rescaled to the total 
branching ratio of the states in (F, R) : 



(F,R) 



n 



(28) 



The scarcity of precise data makes it difficult to disentangle the various SU(3) contributions. 
A very rough estimate of the relevant branching ratios has been made (11) in order to estimate 
y in this approach. The states considered are PP, PV, VV, 3P (three pseudoscalars) , and 4P. 
In most cases, the contributions to y are O(10~ 3 ) at most. The only exception to this is 
the 4P intermediate states, where very large values of y^p^R are found. This dramatic effect 
comes mostly from the fact that some members of the SU(3) multiplets that should cancel 
have run out of phase space, since the average mass of the multiplet is not far from threshold 
to begin with. In this case, it is found that y ~ 10~ 2 can be obtained, although with large 
uncertainties. 

In obtaining this estimate of the long-distance contributions to y, it was assumed that 
phase space is the only source of SU(3) breaking. This constitutes a good estimate if addi- 
tional sources of breaking do not cancel the effect of phase space significantly. However, this 
cancellation would naturally occur if the charm quark is heavy enough so that quark-hadron 
duality, in the spirit of the OPE approach discussed above, is valid. Furthermore, in some 
cases SU(3) breaking occurs in the opposite direction to the one dictated by phase space: 
B(D° — > K + K~)/B(D° — ► 7r + 7r~) ~ 2.9. In any case, the effect found in Reference (11) is a 
reminder of the difficulty of estimating reliably the size of the standard- model effect. 

In this examination of the long-distance contributions to mixing, we have seen that 

x<,10~ 3 , y^lO' 2 . (29) 

The possibility of having y at the 1% level raises a serious question as to whether it will 
be possible to extract a smaller value of x from data with enough precision to constrain the 
short-distance physics contributions from the standard model and/or its extensions. 



2.2 CP Violation 

In the standard model, the D system is not as sensitive to CP violation as the K and B 
mesons are. Once again, the small effects predicted in the standard model could leave open 
a window to the observation of new-physics effects. Here — rather than going into specific 
calculations — we discuss some general features of CP violation in D mesons, both in the 
standard model and beyond. 



2.2.1 DIRECT CP VIOLATION 

Direct CP violation requires the presence of both weak and strong relative phases between 
two or more amplitudes contributing to a given final state. In the standard model, relative 
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weak phases can only be obtained in Cabibbo-suppressed decays, for instance, via the inter- 
ference between spectator and penguin amplitudes. In order to estimate the size of the CP 
asymmetries this would generate, we write 

Acp a IHV.V^K,] sini>t p m 

~ A 2 n A 4 sin | £ HT 3 , 

where <5 st is the strong relative phase between the penguin and the spectator amplitudes, P 
and S, and A ~ 1 and 77 are CKM parameters in the Wolfenstein parameterization. Specific 
model calculations (12) for D — ► KK,irir, K* K , three-body modes, etc. yield this order of 
magnitude for the effect. New physics could enter, for instance, through large phases in the 
penguin diagram. This could give asymmetries of the order of 1% or larger. On the other 
hand, Cabibbo-allowed decays do not have two amplitudes with different weak phases, and 
therefore the CP asymmetry is zero in the standard model. Some new-physics scenarios, e.g., 
some left-right-symmetric models, provide extra phases and could give asymmetries as large 
as 1%. 



2.2.2 INDIRECT CP VIOLATION 

When CP violation is negligible, the CP asymmetry is proportional to 

Acp ~ —2 (xcos<5 + y sin 5) sin</> (Tt). (31) 

In the standard model, (j> ~ 2A 2 X 4 r] ^ 10~ 3 . New physics in x could induce a large sin</>. 
Thus, even if x -C y, provided there is a significant relative strong phase between the Kir final 
states, there could be sensitivity to CP-violating phases from new-physics contributions to x. 



2.3 Constraints on Physics Beyond the Standard Model 

Because the predictions for D° — D° mixing in the standard model are quite uncertain, the 
existence of a discovery window between the current experimental limits and the standard- 
model prediction cannot be established accurately. On the other hand, several extensions of 
the standard model predict large enhancements in Amp,. In many cases, even the current 
experimental limit is enough to severely constrain models. The situation is illustrated in 
Figure 2, where standard-model predictions for x and y are collected along with predictions 
from extensions of the standard model for x. Current experimental bounds already exclude 
or severely limit the parameter space of some standard- model extensions. In some of these 
models, theorists have been forced to choose either theoretical mechanisms or parameters that 
avoid FCNC signals in K and B physics. However, this often leads to large D° — D° mixing 
effects. Thus, although a positive observation of x ^ 10 -2 may not be a clear indication 
of new physics, pushing the experimental bound as low as possible has a great impact on 
theoretical model building. We illustrate this point in two paradigmatic cases: weak-scale 
supersymmetry (SUSY) and strong dynamics at the TeV scale. 



2.3.1 THE MINIMAL SUPERSYMMETRIC STANDARD MODEL and Am D 

Weak-scale SUSY is a possible solution to the hierarchy problem. The minimal supersymmet- 
ric standard model (MSSM), the simplest SUSY extension of the standard model, involves a 
doubling of the particle spectrum by putting all standard-model fermions in chiral supermul- 
tiplets and placing the standard-model gauge bosons in vector supermultiplets. Many new 
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parameters are introduced. The soft SUSY-breaking sector generally includes three gaugino 
masses, as well as trilinear scalar interactions and Higgs and sfermion masses. 

In general, sfermion masses are not related to fermion masses. In particular, if we choose 
to rotate the squark fields by the same matrices that diagonalize the quark mass matrices, 
squark mass matrices are not diagonal (14). In this "super-CKM" basis, squark propagators 
can be expanded so that nondiagonal mass terms result in mass insertions that change the 
squark flavor. These mass insertions can be parameterized in a model-independent fashion 
via 

/ sru \ _ (^j)aA' /qo\ 
^ AA ' _ ~Jtf2 ' ( S2 > 

where i ^ j are generation indices, A, A' denote the chirality, (M^) 2 are the off-diagonal 
elements of the up-type squark mass matrix, and Mq represents the average squark mass. 
This effect can be avoided in specific SUSY-breaking scenarios such as gauge mediation or 
anomaly mediation, but is present in most situations, e.g., if SUSY breaking is mediated 
by gravity. In most cases, the main contribution from these soft SUSY-breaking terms to 
FCNC transitions comes from the gluino-squark loops. Such contributions to rare K and B 
transitions have led to strong universality constraints on the charged Q = — | squark sector 
(see, e.g., (15)). The most stringent bounds that apply to the nonuniversal soft breaking 
terms (S^w come from the experimental searches for D°-D° mixing. 2 For instance, the 
CLEO limit (16) implies (17) 

where 6 is a strong relative phase between the Cabibbo-allowed and the DCS D° — ► Kir decays. 
Neglecting this phase results in the constraints obtained in Reference (17), which we collect 
in Table 1. These bounds were obtained assuming that (Si 2 )rr = and {5\ 2 )lr = {$12) RL- 

Similar bounds on (Sfj)\\i are obtained in K and B processes. In general, SUSY models 
naturally give (5^)ll and (^)ll of order 1, whereas (<5|-)lk — (mz/Mq) is expected. This 
can be avoided if squark masses are assumed to be universal at some high energy scale. 
However, the renormalization group running will generate nondiagonal values of the (%)s in 
general. The mechanism of quark-squark alignment (19), in which squark mass matrices are 
diagonalized along with the quarks, are an alternative to squark degeneracy. This alignment 
occurs naturally in theories with an Abelian horizontal symmetry (20). In order to avoid 
K°—K° mixing bounds, the down-squark mass matrices are chosen to be diagonal. In order to 
obtain the correct Cabibbo angle, an off-diagonal term in the up-squark mass matrices must 
be allowed. This leads to the prediction of large D° — D° mixing, saturating the experimental 
bound. 



2.3.2 CONSTRAINTS ON STRONG DYNAMICS FROM D° - D° MIXING 

We briefly address the constraints arising from Amc on theories where the electroweak sym- 
metry is broken by new strong interactions at the TeV scale (for a review, see (21)). These in- 
clude technicolor (22) and extended technicolor (ETC) (23), top condensation models (24, 25), 
and composite Higgs models (26), among others. In fact, this may also cover the standard 
model if the Higgs is not light enough (27). 

2 Limits obtained from charge and color breaking and bounding the potential from below (18) apply to the 
trilinear terms but not to the squark mass terms. Thus, unless the squark mass matrices are kept diagonal, 
such arguments cannot be used to constrain the nonuniversal mass insertions. 



D°-D° MIXING AND RARE CHARM DECAYS 



11 



If we consider the standard model as an effective theory, then its scalar sector depends on 
two quantities (28): A, the scale of the underlying physics (or scalar "compositeness" scale), 
and /, the amplitude for producing the scalars from the vacuum. The ratio of these two scales 
k = A// is estimated to be k ~ 4tt in naive dimensional analysis (29). On the other hand, 
in order for fermions to acquire masses from their Yukawa couplings to the Higgs, they may 
interact with its constituents through additional interactions characterized by the mass scale 
M. Assuming that this flavor scale is associated with new gauge interactions, it is possible to 
make explicit the dependence of the Yukawa couplings on A and M (27): 

2 A 2 

-TT2 — fR<WL ~ ^—1rH>u (34) 

where tpL is an SU(2)^ doublet containing /j, and v — 246 GeV is the vacuum expectation 
value of the Higgs field. Thus, in order to have the correct fermion mass, one has 



A ^ — .i/\/2rt— . (35) 
9 V v 

The relation in Equation 35 between the Higgs compositeness scale and the scale of flavor 
physics will result in bounds on A. In particular, in addition to the interactions coupling 
the constituents of the Higgs to the standard-model fermions, there are diagonal interactions 
involving only the standard-model fermions leading to dimension-six operators of the form 

g 2 - 

J^2 {fL,RlnfL,R) (fL,RJiMfL,R)- (36) 

For instance, for quarks, the interactions in Equation 36 result in FCNC once the quark mass 
matrix is diagonalized. Thus, operators such as 

■£2 U ™*R Ulr (ul,rWl,r) (u l ,r^c l , r ) , (37) 

where the matrices Ul,r rotate the left-handed and right-handed up-quarks to their mass 
eigenbasis, will induce D° — D° mixing. Similarly, there will be FCNC operators acting in 
the K and B sectors. Bounds on M/g then result in bounds on A. For instance, assuming 
Uf^ R ~ sin#c, one obtains 

A;s22Tev fmF) <38) 

from the LR product of color-octet currents. For the naive-dimensional-analysis value k ~ 4tt, 
we have A £ 78 TeV. 

In principle, this bound could be weakened considerably by assuming Ujj r <C 1. But to 
do so would invite bounds from the down-quark sector. For instance, if Df s ~ sin #(7, then 
A ^ 15 TeV. This is a general result that applies not only to ETC theories but also to a host 
of other scenarios. 

The constraints of FCNCs apply not only to the flavor-physics scale but also to the Higgs- 
compositeness scale. If we interpret the FCNC constraint as a bound on the cutoff scale, 
these theories are fine-tuned. Whatever the way out of the bounds, this analysis suggests 
that these scenarios tend to have plenty of flavor physics. Furthermore, we can interpret the 
cutoff scale A as an implied upper bound on the Higgs mass by making use of the triviality 
constraint (27): 

2 , / A \ 4vrV , . 

m » h, (™) s — • (39) 
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The typical flavor-physics bounds on A ^ 80 TeV yield run < 400 GeV. Then, if we consider 
the standard model as an effective theory valid up to a scale A, a light Higgs results in a 
more fine-tuned scenario, i.e., m#/A <C 1. In heavier-Higgs scenarios, the cutoff scale is lower 
and flavor physics is likely to saturate the experimental bounds on a variety of flavor-physics 
observables. This circumstance is typical in theories where the Higgs is a composite field and 
the standard-model fermions feel the interaction binding it. On the other hand, these bounds 
can be avoided in models with a GIM mechanism built in (30). 

Topcolor: In top-condensation models (24), the constituents of the Higgs are the third- 
generation left-handed quarks as well as tR. Hill (25) proposed that a new gauge interaction 
strongly coupled to the third-generation quarks is responsible for top condensation. The 
topcolor interactions break at the TeV scale as SU(3)iX SU(3)2 — > SU(3) co i or , leaving, besides 
the massless gluons, a set of color-octet gauge bosons (the top-gluons) leading to the Nambu- 
Jona-Lasinio effective interactions that result in top condensation. This leads to electroweak 
symmetry breaking as well as to a large "constituent" top mass. 

Tilting the vacuum in the top directions to avoid a large 6-quark mass is typically accom- 
plished through additional Abelian interactions that leave a Z' strongly coupled to third- 
generation fermions. In some models, the tilting is done by simply arranging that bn not 
couple to the topcolor interactions. The top-gluon interactions (as well as the Z's if present) 
are nonuniversal, leading to FCNC at tree level. These arise after quarks are rotated to their 
mass eigenbasis by the rotations 

where the rotation matrices Ul,r and are unitary. The CKM matrix is then Vckm = 

u\v L . 

Constraints on topcolor models are reviewed in Reference (31). The bounds from the down- 
quark sector impose severe constraints on the entries of ~Dl,r mainly coming from the exchange 
of bound states that couple strongly to the b quark. There are several contributions to Amp. 
First, if we consider theories where topcolor only partly breaks the electroweak symmetry, 
there will be a set of pseudo-Nambu-Goldstone bosons in the spectrum with masses of a few 
hundred GeV. They contribute (31) 

^^(^^urutf, (4i) 

where f% ~ (50-70) GeV is the top-pion decay constant. If it is assumed that the entries of 
Ul,r are similar to those of Vckm, we then obtain 

A TO ^xH-" f 2 °° GeV ) 2 GeV, (42) 

which is just below the current experimental bound even if we consider 5 = 0. 

However, the scalar sector of topcolor theories is somewhat model-dependent. For instance, 
in some scenarios there are no top-pions in the spectrum and topcolor interactions fully break 
the electroweak symmetry. The most model-independent aspect of topcolor theories is the 
top-gluon sector. The constraint that the third-generation coupling be supercritical leaves the 
top-gluon mass as the only free parameter, beyond the quark-rotation matrices. Top-gluon 
exchange also contributes to B°-B° (32) mixing, from which important bounds on topcolor 
models can be derived. 

Here we examine the top-gluon contributions to D° — D° mixing. In principle, many four- 
quark operators, upon quark rotation to the mass basis, contribute to Amp as a result of 
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top-gluon exchange. Before rotation, we have the product of third-generation currents, which 
is enhanced by a factor of (31) cot 2 9 ~ 22; the product of one third-generation current and 
one light current, proportional to cot 9 x tan# = 1; and the product of two light currents, 
which is suppressed by tan 2 9. It appears at first that — although suppressed by tan 2 9 — the 
last group would result in the largest contribution, since the quark rotations in the first two 
generations are expected to be of the order of the Cabibbo angle. A typical contribution to 
Amn from the last group gives (99) 

Am D ~ 4ira s (M G ) tan 2 0^^~ (u L 7^u L )(u L ^u L ), (43) 

G 

which by itself gives a very stringent constraint, Mq/U^ 1 > 100 TeV. 

This inference is not necessarily correct. Not only are the matrices Ul,r unitary, but it is 
also very plausible that the 2x2 block of Ul,r rotating the first two generations is nearly 
unitary. This would be the most likely situation if both Ul,r have a form similar to Vqkm- If 
this is the case, then 

7 JUU*1 JUC 7 /CU*l ICC { A A\ 

u l,rUl,r - -Ul,rUl,r- ( 44 ) 

This implies that all contributions involving one or two light currents before rotation will al- 
most cancel, thus leading to no real bound. We conclude that the only significant contributions 
to D° — D° mixing come from 



Am D ~ 4™ S (M G ) ^ {U^rU^rY (u l ^u l )(u l ^u l ). (45) 

G 

2 



Using f D = 0.2 GeV, a s (l TeV) = 0.09, and cot 2 9 = 22, this leads to 

„ {Uf u * R Uf c R \ 2 (\ TeV\ 2 
Am D ~ 1.8 x 10- 13 L :\ ^ R ^j- GeV, (46) 



sin 5 # c J \ M G J 

where we assume that the elements o£Ul,r scale with the Cabibbo angle Qq ~ 0.22 just as the 
elements of Vckm do. This is compatible with the current experimental limit for Amo, e.g., 
for Mq > 2 TeV or with some of the elements of Ul,r being slightly smaller that the naive 
scaling. Thus we see that generically the contributions of top-gluons, although not in conflict 
with bounds from D° — D° mixing, are not far from the current experimental sensitivity. 



2.3.3 OTHER NEW-PHYSICS SCENARIOS 

As is obvious from Figure 2, there is no shortage of theories beyond the standard model 
whose parameter space is either already constrained by the current experimental limit on 
D° — D° mixing or sits just below it. Here we mention just a few more. 



Extended Higgs sectors without natural flavor conservation: In these scenarios, 
scalars have tree-level FCNCs that may saturate the current limits on mixing just as in the 
case of topcolor scalars. 



Fourth-generation down quarks: Whether a fourth-generation down quark is a member 
of a doublet or a singlet, its contributions to the box diagram dramatically break the GIM 
cancellation in Equation 10 if it has a large Yukawa coupling to the Higgs field. It could 
saturate the bound if its mixing with the standard- model quarks is judiciously chosen. 
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Extra dimensions: If gauge or matter fields propagate in compact extra dimensions, the 
low-energy effective theory below the compactification scale 1/R will contain the massless 
states, the zero-modes. Starting at the compactification scale there are massive Kaluza-Klein 
(KK) excitations of these fields. These towers may contribute to flavor observables and in 
particular to D — D° mixing through loops. In addition, if fermions propagating in the bulk 
have different bulk mass terms, their wave function in the extra dimensions will differ (see, 
e.g., (33)). This leads to nonuniversal couplings to the KK excitations of the gauge fields 
and therefore to tree-level FCNC interactions. For instance, the interactions with the first 
KK excitations of gluons may mimic those of top-gluons if the third-generation quarks are 
localized toward the TeV brane in a compact five-dimensional model with Randall-Sundrum 
metric. Of course, as discussed above, the observation of any nonvanishing value of x cannot 
be interpreted as evidence of new physics because of the proximity of the highest estimates 
from long-distance, hadronic physics. However, this observation may be complementary to 
any deviation from the standard model observed in other flavor-physics observables as well as 
at the energy frontier. 

2.4 Experimental Status and Prospects 

The charm quark was discovered in 1974 in e + e~ annihilation at SPEAR at the Stanford 
Linear Accelerator Center and simultaneously in hadronic collisions at Brookhaven National 
Laboratory (34). Since then, charm hadrons have been produced by a wide range of exper- 
imental facilities. The principal production processes are e + e~ — ► cc at y/s = m[T(45)]; 
Z° — ► cc; hadroproduction, both at fixed-target experiments and at the Fermilab Tevatron; 
photoproduction; and threshold production, e + e~ — ► ^(3770) — > DD. The cross sections vary 
from 1.3 nb in e + e~ — > cc at y/s = m[Y(4«S)] to microbarns for photoproduction, and to of or- 
der a millibarn at the Tevatron. However the ratio of signal cross section to background cross 
section varies from 1/2.5 in e + e~ annihilation at the B factories to 1/100 at the Tevatron. 

Trigger efficiencies vary by orders of magnitude between experiments, and detection effi- 
ciencies vary as well. In a typical (30 fb _1 ) year at a B factory, 8 x 10 7 charm quarks are 
directly produced, with a similar number produced through the decay of B mesons. The 
trigger efficiency is close to unity, and the detection efficiency is of order ~ 50% for a two- 
body hadronic decay. Because of the high relative production rate of charm compared to 
background, excellent mass resolution, and excellent particle identification, there is no need 
to exploit the long lifetime of the charm hadrons to isolate a clean charm data sample. At 
the Z°, the environment is similar but the charm hadrons have higher momenta, ~ 30 GeV/c 
compared to 3 GeV/c near B threshold (35). In contrast, fixed-target experiments measure 
the c-hadron decay time very precisely, which is crucial to isolate clean event samples from the 
large noncharm backgrounds. The early part of Run Ha at the Tevatron produced ~3x 10 11 
charm quarks. Until recently, most of these particles would not have been observed. However, 
there was a tremendous breakthrough in 2002, when CDF (36) demonstrated the ability to 
trigger on the detached vertices produced by hadronic B- and D-meson decays (36). Some 
450,000 D° -> K~ir+ decays have now been reconstructed. 

To give another idea of how significant an achievement this is, it is useful to define an 
effective reconstruction cross section for a given mode at a given experiment as the quantity 
which, when multiplied by the integrated luminosity, yields the number of reconstructed 
events in that mode for that experiment. For CDF the effective reconstruction cross section 
for D° — ► K~ir + is about 7 nb. This should be compared to the effective reconstruction cross 
section of about 13 pb for D° -> K~ir + in e+e" -> D°X at 10 GeV. 

The Cabibbo-allowed two-body decay D° — > K~ir + is easy to reconstruct experimentally 
and has a large branching ratio, 3.85% (9). The number of reconstructed D° — > K~ir + decays 
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therefore serves as a benchmark to quantify the size of charm data samples and facilitates 
comparison between experiments. Using this benchmark, Table 2 compares experiments 
that have recently published results on charm mixing, rare decays, and searches for CP 
violation. Table 2 also lists the time resolution for these experiments, which is important 
in the measurement of lifetimes and mixing. Fixed-target experiments have excellent vertex 
and proper time resolution, five to eight times better than those of the e + e~ experiments 
operating near B-B threshold (for a discussion of the measurement of proper time in heavy- 
flavor physics, see (43)). At the Tevatron, the proper time resolution is similar to that achieved 
in fixed target experiments. 

These data samples have led to many beautiful measurements of the properties of charm 
hadrons (42). The lifetime of charm mesons is now known with exquisite precision; for ex- 
ample, the lifetime of the D° and D + are known to 3 and 6 per mille, respectively. However, 
the data samples collected so far have been too small to allow detection of charm mixing, 
CP violation, and (with one recent exception) rare decays. Although the fixed-target pro- 
gram at Fermilab is now complete, the outlook for much larger samples of charm hadrons 
in the near future is very promising. The experiments BaBar, Belle, and CDF expect to 
increase their data samples, compared to those used in current publications, by factors of 
5, 10, and 30, respectively, over the next few years. Each of these experiments also has 
an upgrade path. There are plans to increase the luminosity of the KEK-B accelerator to 
~ 10 35 cm _2 s _1 by about 2007, an order of magnitude greater than at present. The config- 
uration is called Super-KEK-B. (The plans and current status of this project are described 
at http://www-kekb.kek.jp/SuperKEKB/Workshop.html.) This impressive luminosity corre- 
sponds to a production rate of 2 x 10 9 charm hadrons in a Snowmass year of 10 7 s. A five -year 
program is envisaged at KEK-B. The SLAC group is developing a proposal for a "super B 
factory" with a luminosity goal of ~ 10 36 cm _2 s _1 , corresponding to 2 x 10 10 charm hadrons 
m 10 7 s. This requires a new machine and a very significant upgrade of the BaBar detector, 
known as SuperBaBar (see, e.g., (44)). 

At CDF there is great potential, if charm stays within the trigger bandwidth, to reconstruct 
30 million D — > K~7r + in the 4.4 fb~ l expected in Run II (our estimate is based on the first 65 
pb~l of integrated luminosity), which would be the largest charm sample so far recorded. The 
sample could grow by another factor of five during Run lib. The success of the CDF detached 
vertex trigger augurs well for future hadron collider experiments such as the proposed BTeV 
experiment (see, e.g., (45, 46)), which expects to collect ~ 8 x 10 s charm hadrons in 10 7 s. 
At the Large Hadron Collider (LHC), the charm cross section is twice as large as at the 
Tevatron, but most of the increase occurs at inaccessibly small angles to the beam direction. 
Detector background rates also increase. The design of the trigger for the dedicated S-physics 
experiment LHCb (see, e.g., (47)) is not expected to have significant acceptance for charm. 
At the LHC experiments, ATLAS and CMS, the emphasis is placed on triggers for high pt 
physics. The experiments do not have detached vertex triggers, although they may be added 
as upgrades a few years after turn on. No projections for charm yields exist at this time (48). 

The prospect of a multi-hundred-GeV e + e~ linear collider has led to an investigation of the 
possibility of running at the Z pole, where the charm-quark cross section is ~ 6 nb. Although 
the linear collider would run well above the Z pole, a scheme exists to allow continuous Z 
running simultaneously, a so called Giga-Z machine (see, e.g., (49)). The production statistics 
of Giga-Z do not, in our view, make it competitive with other future facilities, so we do not 
discuss it here. 

The reaction ^(3770) — > DD occupies a special place in the charm experimentalist's and the- 
orist's arsenal. It is the only charm production process in which the charm-quark/anticharm- 
quark pair produces two charm hadrons that are quantum-mechanically correlated. It is also 
the only experimental environment in which the probability to reconstruct two charm hadrons 
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in a single event is large. The "0(3770) therefore offers crucial experimental advantages for the 
determination of absolute charm branching ratios, charm-decay constants and charm semilep- 
tonic form factors, and the CKM matrix elements V cs and V C( i- This suite of measurements 
is important to the international program in precision flavor physics (see, e.g., (49)) and is 
widely held to be the main motivation for a charm factory. Particularly pertinent to this re- 
view, the "0(3770) offers unique opportunities to search for charm mixing and CP violation by 
exploiting quantum coherence and to search for rare decays by exploiting a background-free 
environment. 

Very limited data have been taken at the "0(3770) since 1984, when Mark III (50) ran there 
and accumulated 9.6 pb~l. More recently, the BESII experiment obtained 20 pb~l at and 
nearby the resonance (51). At the "0(3770), CLEO-c (52) expects to accumulate a data sample 
300 times larger than MARK III, whereas BES III 3 will have a sample 10 times larger than 
CLEO-c. To demonstrate the power of 0(3770) running, we also consider the reach of a purely 
hypothetical "Super Charm factory" with a luminosity 17 times larger than that of BEPC II, 
although it is considered unlikely that such a machine could be built. Table 3 summarizes the 
future charm data samples we use for projections in this review. 

In the kaon and B® meson systems, the mixing rates are large, (Am/T) ~ 0(1). The mixing 
parameters can be measured directly from observation of flavor oscillations as a function of 
time. In the B® system, the mixing rate is expected to be so large in the standard model that 
excellent time resolution is required to detect the oscillations. The situation is different for D 
mesons, where the mixing rate is expected to be very small, precluding the direct observation 
of flavor oscillations. Alternative methods include the following: 

1. Measurements of wrong-sign D° decays, using either semileptonic final states or hadronic 
final states. 

2. Comparison of the lifetime of the D° measured in hadronic decays and to final states 
that are CP eigenstates. 

3. Time-integrated measurements exploiting quantum coherence at the 0(3770). 

We now consider each of these methods in turn. Only the first two have been attempted so 
far. A detailed discussion of recent results can be found elsewhere (53). 

In order to detect mixing, which is a small effect, it is necessary to know the initial flavor 
of the D meson and the flavor at the time it decayed. In most measurements that have been 
performed so far, both at fixed-target and at e + e~ facilities, the initial flavor of the D° is 
determined from the observation of the sign of the "slow pion" (also called the "soft pion") 
in the decays D* + — ► D°irf and D*~ — > D°ir~, where s denotes soft. This strong decay 
is extremely useful in charm and beauty experiments because the branching ratio is large 
D* + — > D°7T S = 68.1%, and Q, the mass difference between the initial and final states, is 
small— Q = m(D* + ) - [m(D°) + m(ir s )] = 6 MeV. Consequently, D* + -> D°tt + produces 
a narrow peak in the Q distribution. In addition, the background is extremely low, so the 
miss-tag rate is small, about one per thousand in current experiments. 

2.4.1 SEMILEPTONIC FINAL STATES 

In semileptonic decays, Af = Aj = and we can write (in the small mixing limit) 



r(t) 



P 



P -Tt 

\9-{t)\ 2 



P 



2 (x 2 + y 2 )(Tt) 2 . (47) 



The plans and current status of this project are described at http://bes.ihep.ac.cn/besIII/index.html. 
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The value of f(t) is obtained by exchanging q and p in Equation 47. In the absence of 
CP violation, the integrated mixing rate is Rm — {x 2 + y 2 )/2. The charge of the soft 
pion tags the flavor of the initial state as a D° or D°. The Cabibbo- favored decay [with 
B(D° — > K~l + v) ~ 7%] gives a (pif£~) "right-sign" (RS) charge correlation. If the D° mixes, 
before it decays it is a source of "wrong-sign" (WS) decays: D° — > l) — > K + l~u. In the 
standard model, only one diagram exists for semileptonic decays and so there is no direct decay 
route for D° — ► K + l~v. Accordingly, we describe the measurement as theoretically clean: the 
WS pair (ir^£~) is unambiguously a mixing signature. However, the time distribution of the 
mixing signal is quadratic in the mixing variables. Therefore, although the (Ti) 2 e -r * term 
suits fixed-target experiments, the measurement is relatively insensitive to small values of x 
and y. Furthermore, x and y are not measured separately. Finally, the undetected neutrino 
complicates the measurement. 

At fixed-target experiments, beginning with E791, the technique is to reconstruct the 
semileptonic decay using m(ir s K + £~) — m(K + £~). The reconstructed peak has a large width 
due to the missing neutrino. However, the detached vertex leads to a very clean signal in both 
e and fi modes. The E791 result is (54) R M = \{x 2 + y 2 ) < 0.5% at 90% CL. CLEO obtains 
Rm < 0.87% @ 90% CL) (55) measured in the channels D° K+fj,'^ and D° -> K* + e~u e , 
respectively. The most restrictive measurement of Rm comes from FOCUS (56) in the mode 
D° K+n-Up. 

Rm ~ \{x 2 + y 2 ) < 0.12% at 90% CL (48) 

The bound on Rm is a circle in the x — y plane centered at the origin. In Figure 3 this bound 
is shown ELS cL circle with horizontal shading. 

The future outlook for searcher of charm mixing using semileptonic final states is as follows. 
B factory results are expected soon. The SuperBaBar working group estimates that, with a 
10 ab _1 data sample, the sensitivity would be Rm < 5 x 10~ 4 (57). The technique can, in 
principle, also be employed by future experiments at hadron machines because the presence 
of the lepton in the final state helps triggering. However, the technique has the greatest 
sensitivity at e + e~ machines operating at the ^(3770). In all cases, if mixing is detected in 
this channel and y turns out to be larger than or comparable to x, a separate measurement 
of y will be needed. 



2.4.2 HADRONIC FINAL STATES 



A D° can produce a WS hadronic final state either by undergoing a doubly Cabibbo-suppressed 
(DCS) decay or by first oscillating into a D° that subsequently undergoes a Cabbibo-favored 
(CF) decay. The WS decay includes three components: one from the DCS decay, a second 
from mixing, and a third from the interference between the first two. Assuming CP conserva- 
tion and expanding the decay rate up to 0(x 2 ) and 0(y 2 ) results in the following expression 
for the time evolution of the hadronic WS decay rate: 



r{t) 



-rt 



Rd + y/R~Ey' (rt) + 



(x> 2 + y' 2 ) 



(Tt) s 



(49) 



whereas the RS decays have simple exponential time dependence oc exp r *. Here Rd is defined 
by 

= -^R- e ~ lS (50) 

and gives the rate for the DCS component. The mixing contribution is quadratic, and hence 
very small, but the interference term is linear in y' and may result in a measurable deviation 
from a pure exponential decay. 
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The specific case / = K + tt has been studied. The parameters x' and y' are related to the 
mixing parameters x and y by a rotation (58) 

x' = xcosSkk + ysinSicK, y = y cos Skw — x sin 8k<k- (51) 

A strong-phase Sktt exists between the DCS and CF decay amplitudes (59, 60). Although 
it vanishes in the SU(3) limit, this can be a badly broken symmetry. In a measurement 
based on the time evolution of the WS state, it is not possible to determine the phase. Since 
the expression is quadratic in x' , its sign is also not determined. We discuss techniques to 
determine 5ktt in Section 2.4.5. 

To search for CP violation in the time evolution of the WS hadronic state one applies 
Equation 49 to D° and D° separately. One determines {-R^S' x %i V+} ^ or ^° candidates and 
{R ws , x'^,y'_} for D° candidates, where R^s are the corresponding wrong-sign rates. The 
separate D° and D° results can be combined to form the quantities 



{ ° ~ r+ i p- ' Am ~ p+ _j_ p- 
a P n P "■M M 



2—^- A M = „T (52) 



where R M = (a;±-|-y±)/2. Here Am parameterizes CP violation in the mixing amplitude and 
leads to 

^(l + ^e**, (53) 
P 

where 4> is responsible for CP violation in the interference between the DCS decay and mixing. 
Also, Ad is related to CP violation in the DCS decay amplitude and is defined by making the 
replacements \[Rp — ► y/Rp(l + Ajj) for the case of initial D° and ^/R~B -> v/flF/(l + Ar>) for 
the D°. In this way, the measurement of the time dependence is sensitive to the combinations 
y' cos <fi± x' sin through 



x± 



'1±A M 



1 T A 



M 



(x' cos cf) ± y' sin (ft) , (54) 



y± = \7 t~t^ (#' cos ( l )±x ' sin ^) • ( 55 ) 

V 1 T ^.M 

An offset in <p of ±7r can be absorbed by a change in sign of both x' and y', effectively swapping 
the definition of the two physical D° states without any other observable consequence. In 
order to avoid this ambiguity, we use the convention that \4>\ < n/2. CLEO uses a different 
approximation (61) in its analysis, which is valid for Ad, Am <C 1 : R% = (1 ± Ad)Rd, £± = 
\/T± Am(x' cos 4>±y' sin <f) , y'± = \/T± Am(j// cos x' sin 0) . 

The total time-integrated hadronic WS rate, assuming CP conservation and normalizing 
to the total RS rate, is 



^ws = jf^y dt =Ro + + 1/2( + y )- (56) 

In this approximation, the mixing rate is Rm = l/2(a/ 2 + y' 2 ) = \/2{x 2 + y 2 ). If there is no 
mixing in the D° system, Equation 56 reduces to Rws = Rp- 

CLEO (61) and BaBar (62) have obtained results with 9.0 fb _1 and 57.1 fb _1 of data, 
respectively. A result from Belle should soon be available (63). All three experiments have 
results for the time- integrated WS decay rate Rws (61, 62, 64). 

WS candidate events of the types D° — » K + ir~ and D° — ► K~tt + are selected by requiring 
the ir s from the D* decay and the daughter K of the D° to have identical charge (WS tag) . To 
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detect a deviation from an exponential time distribution in WS events, a likelihood fit to the 
distribution of the reconstructed proper decay time t is performed. The likelihood fit includes 
a signal and a background component and models each as the convolution of a decay-time 
distribution and a resolution function. 

The WS sample is statistically limited, amounting to ~ 1/300 the size of the RS sample 
selected with the same criteria except for the WS tag. The RS sample is used to constrain 
aspects in the fit that are common to the two samples. CLEO and BaBar determine the 
resolution functions of WS and RS signals and of the common background types with the RS 
sample. In the WS sample, a significant additional complication arises from the much lower 
achievable purity, e.g., ~ 50% in the CLEO analysis. 

CLEO and BaBar both fit using Equation 56 and consider the case with and without CP 
violation. The results of both experiments are consistent with the absence of mixing and 
CP violation. We first consider the case with CP conservation. CLEO finds one-dimensional 
limits of (l/2)x' 2 < 0.038% and -5.2% < y' < 0.2% at the 95% CL. The systematic error is 
±0.2% (±0.3%) for x' [y') dominated by knowledge of the background shapes and acceptances. 

CLEO's two-dimensional result is a contour on the x'—y' plane that contains the true value 
of x' and y' with 95% confidence. The contour is constructed using a Bayesian approach. 
Systematic errors are small but are not included. The contour is shown as the cross hatched 
region in Figure 3, where the strong phase shift Sk-k between the Cabibbo-favored and DCS 
decays has been assumed to be zero so that the range in x' (y') corresponds to the range of 
the allowed region in x (y). 

In their fit, the BaBar collaboration allows x' 2 to take nonphysical negative values and 
finds that the most likely fit point when CP is conserved has a negative value of x' 2 . BaBar 
observe a correlation between x' 2 and y' . They present their results as a contour on the x' 2 —y' 
plane (see (64)), where the contour has been constructed using a frequentist approach based 
on toy Monte Carlo and systematic errors have been included. BaBar made their contours 
available to the authors of this review and we have redrawn the x' 2 -y contour on the x-y 
plane by converting x' 2 to xi and assuming 5k-k to be zero. The BaBar redrawn contour is 
displayed in Figure 3. The numerical result is given in Table 4, where the BaBar limits on 
x' 2 and y' are obtained by projecting their contour onto the corresponding axis. Also shown 
in Figure 3, and numerically in Table 4, is a preliminary contour from FOCUS (65) derived 
from WS D° -> K+-k~. 

The BaBar upper limit for x' 2 is almost three times larger than CLEO's but is based on a 
data sample six times larger. However, CLEO and BaBar have used different techniques to 
obtain the 95% CL limits, and the treatment of the fit output parameter x' differs as well, 
as can be seen in Figure 3. The region allowed by the BaBar result is about a factor of two 
more restrictive than that from CLEO. This is to be expected, since the measurement of x (y) 
using this technique scales as 1/L 1 / 4 . For a discussion with an alternative point of view see 
Reference (53). The sensitivity of this technique to x and y will improve by a factor of two for 
each B factory with 500 fb . 

A direct comparison of the CLEO and BaBar results is not possible when CP violation is 
allowed in the fit, because CLEO uses as fit output parameters x', y', Rd and Ad, Am, and 
sin0, whereas BaBar uses Ad,x' 2 , y' + , R^, R M (x f2 , and y'_, R D , R M ) for the D° (D°) case. 
Fits that allow for CP violation lead to slightly less restrictive limits on x' and y' (se Table 4). 
In no case is evidence for CP violation found, i.e. BaBar finds Ad and Rm consistent with 
zero, and CLEO finds Am, Ad and sin</> consistent with zero. 

Extraction of the time-integrated WS decay rate i?ws : CLEO and BaBar measure 
-^ws by repeating the fits described in the previous section with the assumption of no mixing in 
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the D° system, i.e., x = y = (61, 62). CP violation is allowed, so the fit returns Rws — Pd 
and Aj^. Belle uses a fit in the Mp — 5m plane to determine the time-integrated number of 
signal events in the candidate samples (64). The ratio of the number of signal events in the 
WS and RS candidate samples yields i?ws> where the systematic error is dominated by the 
uncertainty on the background shapes used in the fit. Table 5 compares results from CLEO, 
BaBar, and Belle to earlier measurements of WS decays by E791 (66), ALEPH (67), and 
FOCUS (68). 

The results are in reasonable agreement. We compute the world average to be 

(R ws ) = (0.368 ± 0.021)%. (57) 

The WS rate is at roughly the level expected in the standard model (see Equation 24). If 
mixing were present, the WS integrated rate would differ from the standard-model expecta- 
tion. 



D° Decays to Multibody Final States: CLEO has also measured Rws m the multibody 
channels D° -> K + tt--k° and D° -> K + tt--k + tt- with the results R ws = (0.43^;^ ± 0.07)% 
and P W S = (0.41 ± 0.04)%, respectively (69, 70). P W S need not be the same for different 
decay modes, but within the errors, i?ws is the same for all decay modes measured so far. 
With the large data samples from the B factories, it may be possible to set D°-D° mixing 
limits using combined Dalitz plot and proper time fits in multibody modes. These modes may 
prove useful in searching for CP violation and understanding strong-phase shifts (71). 

The decay D° — ► K f gir + -K~ may be used to measure x and y directly, since the strong-phase 
difference may be extracted simultaneously in a time-dependent fit to the Dalitz plot. This is 
possible because both the RS and WS decays in the submode D° — ► K*^tt^ have the same 
final state. Thus, one can fit for the phase difference directly. The sign of x can also be 
extracted from such a fit. CLEO presented evidence for a WS amplitude and measured the 
branching fraction relative to the RS mode to be 

|££f^l_ (0.5 ± 0.2«;5 ±i g;?)% (58) 

and the strong-phase difference between the RS and WS to be (189° ± 10 ± 3I5 5 ) using a 
time-independent Dalitz plot fit (72). The last uncertainty is due to the choice of resonances 
and model. No CP- violating effects were observed when the sample was separated into D° 
and D° subsamples. Results of a time-dependent fit with limits on x, y, and CP violation are 
expected soon from CLEO, BaBar, and Belle. This channel may offer the greatest sensitivity 
to x at the large integrated luminosities already collected by Belle and BaBar. 



2.4.3 THE MEASUREMENT OF y 

The mixing parameter y can be determined by measuring the lifetime difference between D° 
decays to CP-even and CP-odd final states. 

T(CP even) - T(CP odd) T(D° K + K~) 
VCP ~ r(CP even) + T(CP odd) ~ T(D° K"tt+) ~ ' (59) 

which then results in 

VCP = V cos 4> - x sin <p (A M + A prod ) , (60) 
where the production asymmetry is defined as 

_ Njpf) - N(D°) 
Apmd - N(D°) + N(D°) ' 1 j 
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It is assumed that (Am, A pro a) <C 1, where Am and <f> are defined as in the previous section. 

Then yep is determined from the slope of the decay-time distributions in samples of D° — ► 
K~~ir + , which is an equal mixture of CP-even and CP-odd final states, and D° — > K~K + or 
7r~7r + , which are even final states. An unbinned maximum likelihood fit to the distribution 
of the reconstructed proper decay time t of the D° candidates at e + e~ machines, and the 
reduced proper time at fixed-target experiments, is performed. Because yep is measured 
from the ratio of lifetimes, many systematic effects cancel. 

A fixed-target experiment, FOCUS, was the first to make this measurement and found a 
value of yep of about 3% that was several standard deviations from zero (68). The e + e~ 
experiments BaBar, Belle, and CLEO have now all made this measurement with greater 
precision than FOCUS. 

The three experiments determine y from unbinned maximum-likelihood fits to the distri- 
bution of the reconstructed proper decay time, t, of the D° candidates. Small biases may 
occur in the measured lifetimes; these are at the 1% level or below. The precision with which 
they are known is limited by the statistics of the simulation and is the dominant source of 
systematic uncertainty for BaBar, whereas for Belle and CLEO knowledge of the background 
contributions to the signal is the largest source of systematic uncertainty. 

The technique, resolution, and systematic errors at y/s = 10 GeV are quite different from 
those in fixed-target experiments. D° candidates are selected by searching for pairs of tracks 
with opposite charge and combined invariant mass near the expected D° mass. The intercep- 
tion point of the D° momentum vector with the envelope of the interaction point (IP) provides 
the production vertex of the D° candidate. The proper decay time of a D° candidate is derived 
from its mass and flight length. 

BaBar and CLEO employ a D* tag. They refit each tt s candidate track with the constraint 
that it coincide with the D° candidate production vertex. This reduces substantially the 
mismeasurement of the tt s , momentum caused by multiple scattering. Then dm = m(D* + ) — 
m(D°) is required to be consistent with the known value. Belle does not use a D* tag but 
requires that the D° candidate flight path be consistent with originating at the IP. All three 
experiments reject events with secondary charm production from P-meson decays with a 
momentum cut. The analyses by BaBar and Belle are for a subsample of the full data sets. 
Table 6 lists the subsample sizes and the number of events. 

Recently BaBar reported an improved analysis that allows for CP violation. If CP violation 
occurs, the D°(t + ) and the P°(r~) will have different lifetimes to decay to CP-even states. 
The effective lifetimes can be combined to form 



y = — -1, (t) = (t + + t-)/2 (62) 
( r ) 

(r) (r++r ) 

Both Y and Ay are zero if mixing is absent. Otherwise, in the limit of CP conservation in 
the mixing amplitude (and little production asymmetry), Y = ycos</> and Ay = xsm<j). All 
BaBar data (91 fb^ 1 ) have been used. Four independent samples were isolated, three tagged 
with a D*: K~tt + measures r°; K~ K + measures (r), and A T ; 7r + 7r~ measures (r), and A T ; 
and an untagged sample of ~ 146,000 K~ K + measures (r). An unbinned likelihood fit to 
mf)0,t,at yields {t),tq, and A T . The statistical uncertainty is small; for example, for K~tt + , 
it is w 0.9 fs, corresponding to 0.5% in y. Assuming the same value of (j), the modes K + K~ 
and 7r + 7r~ can be averaged. BaBar combines the modes to find 

y = (0.8 ± 0.4±£|)% (64) 
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AY = (-0.8 ±0.6 ±0.2). (65) 

This is the single most precise value of Y and the first measurement of Ay. 4 Because Y 
and Ay are consistent with zero, there is no evidence for mixing or CP violation. 

In the limit of CP conservation Y = y. We average Y with previous measurements of 
yep- Figure 4 summarizes the experimental situation: the FOCUS measurement is high, 
whereas the recent BaBar and Belle results move the world average closer to zero. The 
new measurement of Y continues this trend. Within the errors, the measured values are 
consistent. To compute the world average, we exclude the first BaBar measurement because 
the new measurement of Y supersedes it. We find 

(VCP) = (0.8 ±0.5)%. (66) 

The average is consistent with the standard-model expectation of a value of \y\ close to zero. 
Figure 3 shows the world-average value of yep on the x-y plane. 

We note that Ay is not systematically limited and so the precision with which it is deter- 
mined can be expected to improve by a factor of y/h with a 500-fb _1 data sample at BaBar 
or Belle. With the 10 ab -1 sample expected at SuperBaBar, the collaboration estimates the 
following reach: \Y\ < 0.1% and |Ay| < 0.06% (57). 

Figure 3 combines the constraints on the x~y plane for all measurements of mixing from 
hadronic and semileptonic decays discussed so far. The strong-phase difference is assumed to 
be zero for this comparison. It is seen that the BaBar measurement of x' 2 — y' is in reasonable 
agreement with the horizontal band corresponding to the world average of y. If 5 = 40° (the 
estimated maximum in Reference (73)), the elliptical constraints on x and y from D° — > K~ir + 
would be rotated by 5 in a counterclockwise direction, resulting in greater overlap of the CLEO 
and FOCUS measurements with the world average value of y. These results would be even 
more consistent if <5 > 90° (58). The overlap of the BaBar measurement with the y band is not 
significantly improved by rotation. The errors are too large to permit any strong conclusion 
at this time. However, it is clearly very important to obtain a measurement of the strong 
phase. 



2.4.4 MIXING VIA QUANTUM COHERENCE AT THRESHOLD 

At the ^(3770)_ and at the ^(4140), reactions ^(3770) _> D°-D° and ^(4140) -> jD° D° 
produce a D° D° pair in a state that is quantum- mechanically coherent. This enables simple 
new methods to measure the D° — D° mixing parameters (74) in a way similar to that proposed 
in Reference (75). Consider, for instance, D°D° — > (K~ tt + )(K~ tt + ) . The initial D°D° state 
is 

\i>= ^{\D\k 1 ,t 1 )D\k 2 ,t 2 )} + Vc\D°(k 2 ,h)D°(k 1 ,t 1 )}}, (67) 

where rjc = ± is the charge-conjugation eigenvalue of the D°D° pair. At the tp(3770), 
r\C = —1, so the system is in the antisymmetric initial state. The case where the D° decay 
is Cabibbo-favored (CF) and the D° decay is DCS (e.g. D° -> K~Tr + D° K~tt + ) is 
indistinguishable from the the case where the D° decay is CF and the D° decay is DCS. Upon 
adding the amplitudes, the DCS contribution cancels. Therefore, for D°D° pairs produced 
from the decay of the -0(3770), the time-dependent rate is 

T(t) oc ^ e - r (' 1 +^)|^| J B(fci)| 2 | J B(^)| 2 (x±y) 2 (ti -t 2 ) 2 , (68) 

2 9 

4 As this review was going to press, the Belle collaboration announced the preliminary result of a new 
analysis allowing for CP violation in mixing and interference. Using 91 fb _1 of data (82). They find yep = 
(1.15 ±0.69)%. They also measure A r = (-0.2 ±0.6 ±0.3)%, where A r differs from AY by a factor (l + y C p). 
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where B = B K - W + = A(D° — ► K 7r + ), i.e., the decay amplitude. Normalizing to T[D°D° — ► 
(K~ir + )(K + ir~)], which is given by 

r(t) oc le- r ^ +t ^\B K ^ + (h)\ 2 \B K+7T -(k 2 )\ 2 , (69) 
implies for the time-integrated ratio 

(K~tt + )(K~tt + )\ x 2 + y 2 



R 



XK-tt+)(K+it-) J 2 
This is similar to the case of semileptonic final states, 



2|n |2 



We note that the DCS amplitude does not cancel for all decay modes (75). For the case where 
one final state is hadronic and the other semileptonic, we have 

These unambiguous signatures for D°-D° mixing arise because of the quantum coherence of 
the initial state and are analogous to the situation for B°-B° pairs at the T(45). 

The measurement of Rm can be performed unambiguously with the decays ip" — > K~tt + K~tt^ 
and ip" — > ^{KXy^vl^^Xy^v. The hadronic final state cannot be produced from DCS 
decays. This final state is also very appealing experimentally, because it involves a two-body 
decay of both charm mesons, with energetic charged particles in the final state that form an 
overconstrained system. Particle identification is crucial in this measurement because if both 
the kaon and pion are misidentified in one of the two L>-meson decays in the event, it becomes 
impossible to discern whether mixing has occurred. 

The number of RS events in the all-hadronic and hadronic-semileptonic channels combined 
produced in 3 fb _1 is expected to be about 50,000, corresponding to ^/Rm < 1% at 95% CL. 
At threshold, the sensitivity scales as x oc L" 1 / 2 , whereas at a B factory, for the measurement 
of mixing through the time evolution of D° — ► K + ir~, the scaling goes like x oc L -1 / 4 . 
This implies that a 3-fb _1 data sample at charm threshold has similar reach to a S-factory 
data sample of 500 fb -1 . Because the charm-threshold measurement is background-free, it is 
statistically limited rather than systematically limited. At BES III, where the data sample 
is expected to be 10 times greater, the limit will improve to \[R~m <■ 0.3%, but only if the 
particle identification capabilities are adequate. If it were possible to obtain 500 fb _1 at the 
V>(3770), the limit would be JWl < 0.08%. 

The investigation of mixing parameters at the -0(4140) provides valuable complementary 
constraints and a useful consistency check. Note that the final-state (K~tt + ) (K~tt + ) and 
(K~ir + ) from the r\c = +1 initial state have a decay width proportional to y' . Because 
Rm is bigger than x or y, the determination of y' from tjc = +1 initial states is comparable 
in sensitivity to measurements from the higher-statistics rjc = — 1 initial state. 



2.4.5 STRONG PHASES AT CHARM THRESHOLD 

We can also take advantage of the coherence of the D mesons produced at the ip" to extract the 
strong-phase difference 5 between the direct and DCS amplitudes that appears in the time- 
dependent mixing measurements. Because the CP properties of the final states produced 
in the decay of the ip(3770) are anticorrelated, one D state decaying into a final state with 
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definite CP properties immediately identifies or tags the CP properties of the state "on the 
other side." If one state decays into, for example, tt°Ks with CP = —1, the other state is 
"CP-tagged" as being in the CP = +1 state. This allows measurement of the branching ratio 
B{Dcp — > K~tt + ) and cos 5. A triangle relation follows from the definition of Dqp- 

V2A(D CP K~ir + ) = A(D° K~ir + ) ± A(D° K~tt + ). (73) 

This implies 

, B(D± -> K-7T+) , . 

where Rp> is the ratio of the DCS decay to the CF mode. We have used the fact that 
R-D ^ ^/Rd and neglected CP violation in mixing, which could undermine the CP-tagging 
procedure by splitting the CP-tagged state on one side into a linear combination of CP- 
even and CP-odd states, thus requiring time-dependent studies. Both effects, however, are 
negligible. Now, if decays of both D + and P_ are measured, cos<5 can be obtained from the 
asymmetry, 

„ B(D+ -» K~tt+) - B(D~ -► K~tt + ) , P , . 

The asymmetry A is expected to be small. Thus we have 

A(cos(5) ss — (76) 



where iV is the total number of CP-tagged K - 7r + from the C = — 1 initial state. At CLEO-c, 
the total number iV is expected to be about 32,000 in one year run of CLEO-c, leading to an 
expected accuracy of about ±0.05 in cos 6. 

Reference (83) outlines an alternative method to extract 5k-k from a measurement of the 
rates of the DCS and CF decays of the type D — > Ktt. This method requires the determination 
of rate asymmetries in D meson decays to K^ir and Kgir, and makes the assumption that 
the two A7 = 1/2 DCS amplitudes have equal phase. Belle (84) finds that the relevant 
measurements of Kj_, and Kg mesons are possible with a statistical precision sufficient to 
constrain Sk^, but their preliminary measurement is not yet sensitive enough to do so: 

As the data sample grows, this measurement will become interesting, but it will not be as 
precise as the determination of 5k-k at threshold. Because the two measurements have entirely 
different systematic errors, the P-factory measurement will serve as a useful cross check. 

2.4.6 SEARCHES FOR CP VIOLATION IN THE P-MESON SYSTEM 

Searches for CP violation in charm decays have reached accuracies of several percent (85, 86). 
In order to measure an asymmetry, one must know the flavor of the decaying meson. Charged 
D mesons are self-tagging; the neutral D meson requires the above-mentioned D* tag to 
determine whether the decaying particle was a D° or D°. 

The fixed target experiments E791 and FOCUS have made some of the most sensitive 
searches for direct CP violation in the charm system. A subtlety in the measurement is that 
there is a production asymmetry between D + and D~ and between D° and D°, and so these 
experiments normalize all asymmetries to some known CF mode where there can be no CP 
violation — for example, 

Acp(KKir) = ri{D+) + ri{D _ y (78) 
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where 

v(D±) = N{D± - KTK±n±) (79) 
V[ ' N{D±^K^±k±)' 1 j 

Signals are isolated according to (a) the significance of the longitudinal separation between the 
decay (secondary) vertex and production (primary) vertex (L/cr), where L is the longitudinal 
separation and a is the calculated resolution in L; (b) the confidence level of the decay vertex 
fit; (c) the vector sum of the momenta from decay vertex tracks, which is required to point to 
the production vertex in the plane perpendicular to the beam; and (d) particle identification 
and momentum cuts that are specific to a given final state. 

CDF recently presented their first search for CP violation in the charm system using the 
two-track hadronic trigger (88). Using 65 pb _1 of data collected during 2002, about 100,000 
D° -> K-JT+, 8,000 D° -> K~K + and 4,000 D° -> vr+vr" candidates pass selection cuts 
including 5m(D*), the impact parameter of the D°, and the projected decay length of the 
D°. The already impressive statistics yield the most sensitive searches for asymmetries in the 
following modes: 

A CP (D° -► K+K~) = (2.0 ± 1.7 ± 0.6)%, (80) 

A CP {D° -> tt+tT) = (3.0 ± 1.9 ± 0.6)%. (81) 

The dominant systematic error comes from the correction for the charge asymmetry for low- 
momentum tracks in the CDF tracking system. 

CLEO has made numerous searches for CP violation in the charm sector. There is no 
production asymmetry or appreciable detection asymmetry in the T(4S") energy region, but 
statistics are limited. Selected recent asymmetry measurements from all of these experiments 
are tabulated in Tables 7 and 8 and in Figures 5 and 6. We compute the following averages: 

(A C p(D° K+K')) = (0.8 ± 1.2)%, (82) 

(A C p(D° vr+O) = (2-7 ± 1.6)%. (83) 

There is no evidence for CP violation at the current level of sensitivity. 

An entirely different way to search for CP violation is to exploit quantum coherence at the 
tp". The production process 

e + e + ^ ^» ^ D°D° (84) 

produces an eigenstate of CP+, in the first step, since the ip" has J PC equal to 1 . Consider 
the case where both the D° and the D° decay into CP eigenstates. Then the decays 

r - /;/; or /i/i (85) 

are forbidden, where /+ denotes a CP+ eigenstate and /_ denotes a CP— eigenstate. This 
is because 

CP(/4/i) = (-1)' = -1 (86) 

for the I = 1 ip" . Thus, observation of a final state such as (K + K\tt + tt~) constitutes evidence 
of CP violation. Moreover, all pairs of CP eigenstates, where both eigenstates are even or 
both are odd, can be summed over for this measurement. This provides a sensitive way to 
detect CP violation in charm decays and could become sensitive enough to see standard-model 
mechanisms. 

Table 9 estimates the total number of events that would be observed for maximal CP 
violation in a one-year run at CLEO-c. The event samples are not large, but the measurement 
is essentially background free. Moreover, this method is important because it has unique 
sensitivity to the quantum-mechanical phase in the amplitude. This measurement can also 
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be performed at higher energies, where the final state D*°P)*° is produced. When either D* 
decays into a it and a D°, the situation is the same as above. When the decay is D*° — ► 7-D , 
the CP parity is changed by a multiplicative factor of 1, and all decays /+/!_ violate CP (92). 

The coherent nature of the wave function prevents measurements of time-integrated mixing- 
induced CP asymmetries from CP-odd D°D° pairs produced at the ^(3770). Thus, at this 
energy, nonzero time- integrated asymmetries would be a manifestation of direct CP violation. 
The measurement of direct CP violation in D° decays at the ^(3770) uses a different tagging 
technique from that at higher energies because the decay D* + — > D°ir + is kinematically 
forbidden. One D° meson in the event, the tag, is partially or fully reconstructed. One 
example of partial reconstruction is lepton flavor tagging, in which the lepton in the decay 
D° — > £ + X is detected and the charge of the lepton determined. With the flavor of the tag 
established, the flavor of the other D° meson in the event is fixed. Table 10 summarizes the 
expected sensitivity for direct CP violation in a one-year run at the V'(3770) at CLEO-c. 

CP-even initial states, such as those produced at the ^(4140) — > ^yD°D°, are amenable 
to time-integrated CP asymmetries that are nonvanishing and depend only linearly on the 
mixing parameter x. In this case, there is a smaller statistical accuracydue to the smaller 
cross section, and the overall CLEO-c sensitivity is of the order of 3% using lepton flavor tags 
only. Backgrounds are expected to be modest and to have a negligible effect on the measured 
asymmetries, in agreement with previous studies (93, 94). 

The studies described here are only some examples of the many possible search strategies 
(92). Dalitz plot analyses uncover interference effects that are sensitive probes of CP- violating 
phases. Although some of these measurements can be performed with comparable accuracy 
at the B factories, the number of studies that can be performed at threshold is very broad and 
includes some unique measurements exploiting the quantum coherence of the initial state. 

Current measurements of CP violation in D° decay are at the several-percent sensitivity 
level. The limits are considerably worse in D + and D$ decay. There are significant oppor- 
tunities to search for the effects of new physics via the mechanism of CP violation both at 
charm threshold and at the B factories and CDF. 

The most reliable way to compare the reach of future experiments for Aqp would be to 
compile a list of sensitivities estimated by each experiment based on detailed simulations. 
In almost all cases the simulations, or other detailed estimates, do not exist. Instead we 
will use the simplistic benchmark of the number of reconstructed D° — > K~ir + to estimate 
sensitivities. The CP asymmetry scales like 1/ ' \JN(D {) — > K~-k + ). To test this method, we 
note that the experiments E687, E791, and CLEO II, which each reconstructed a few x 10 4 of 
these decays, had a similar CP reach of 5 Aqp ~ 5%. Using this result as a normalization, one 
expects 5Aqp ~ 3% for FOCUS and 5Acp ~ 2% for CDF. Both are in agreement with the 
precision reported by the experiments. We compute the sensitivities for future experiments 
in Table 11. Experiments can be expected to probe CP asymmetries at the 10 -3 level during 
this decade. In the early part of the next decade, we can hope to see CP asymmetries probed 
at the 10~ 4 level. 

3 RARE CHARM DECAYS 

We have seen that charm processes tend to be affected by large nonperturbative effects. 
However, for some modes, a window exists in which theoretical predictions are sufficiently 
under control to allow tests of the short-distance structure of the FCNC transition. This 
is, to some extent, the case in c — > u£ + £~ modes, and therefore we concentrate on their 
potential. On the other hand, radiative charm decays, such as those mediated by c — ► uj, 
are largely dominated by long-distance physics. Their experimental accessibility presents an 
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opportunity to study purely nonperturbative effects in radiative weak decays. We first review 
the standard-model predictions for the leptonic, semileptonic, and radiative decays. Then 
we study the potential for new-physics signals in c — ► u£ + £~. Finally, we survey present 
experimental knowledge and future prospects in rare charm decays. 



3. 1 The Standard- Model Predictions 



The short-distance contributions to the c — > u transitions are induced at one loop in the 
standard model. It is convenient to use an effective description with the W boson and the b 
quark being integrated out as their respective thresholds are reached in the renormalization 
group evolution (96). The effective Hamiltonian is given by (97-99) 



T~t c s 



AG F 



T~t c s 



V2 
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+ Y t Ci<jJ.)O i (j i ) 
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, fj, < m b , 



(87) 



with {Oi} being the complete operator basis, {Cj} the corresponding Wilson coefficients, and 
\x the renormalization scale; the primed quantities are those for which the b quark has been 
eliminated. In Equation 87, the Wilson coefficients contain the dependence on the CKM 
matrix elements V qq i . The CKM structure of these transitions differs drastically from that 
of the analogous i?-meson processes. The operators 0\ and O2 are explicitly split into their 
CKM components, 



0[ q) = {ullA^L^t) 



L> ' 



(88) 



where q = d,s,b, and a, (5 are contracted color indices. The rest of the operator basis is 
defined in the standard way. The QCD penguin operators are given by 



O3 
Or, 



intact) Y,^W R ) » o Q = (ut^4) Ytiskftii)- 

The electromagnetic and chromomagnetic dipole operators are 

7 = j^m c (u L a^c R )F^ , 8 = ^ mc (u L a^T a c R )G^ ; 
and finally the four-fermion operators coupling directly to the charged leptons are 



(89) 
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The matching conditions at jjL = Myy for the Wilson coefficients of the operators 0\-q are 



: (^7/^X^75^) • 



(91) 
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C 3 . 6 (M W ) = 0, C q 2 (M w ) = -A 



(92) 
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with X q = V* q V uq . The corresponding conditions for the coefficients of the operators O 
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(93) 



In Equation 93, we define X{ = mf/My^; the functions -Fi(x), F 2 (x), and C(x) are those 
derived in Reference (100), and the function D(x) was defined in Reference (97). 

To compute the c — > u£ + £~ rate at leading order, operators in addition to O7, O9, and 0\q 
must contribute. Even in the absence of the strong interactions, the insertion of the operators 
0$ in a loop would give a contribution sometimes referred to as leading-order mixing of C 2 
with C9. When the strong interactions are included, further mixing of the four-quark operators 
with O7_io occurs. The effect of these QCD corrections in the renormalization group running 
from M\y down to = m c is particularly important in Cj S (m c ), the coefficient determining 
the c — > wy amplitude. As was shown in Reference (97), the QCD-induced mixing with 2 
dominates Cj S (m c ). The fact that the main contribution to the c — >• wy amplitude comes from 
the insertion of four-quark operators inducing light-quark loops signals the presence of large 
long-distance effects. This was confirmed (97, 98) when these nonperturbative contributions 
were estimated and found to dominate the rate. Therefore, we must take into account effects 
of the strong interactions in Cy ff (m c ). On the other hand, the operator Og mixes with four- 
quark operators even in the absence of QCD corrections (101). Finally, the renormalization- 
group running does not affect 0\q, i.e., Cio(m c ) = C\q(Mw)- Thus, in order to estimate the 
c — > u£ + £~ amplitude, it is a good approximation to consider the QCD effects only where they 
are dominant, namely in Cj S (m c ), whereas we expect these to be less dramatic in C| ff (m c ). 



The leading-order mixing of with O9 results in 



CP cff = C 9 (M W ) + £ Xi 

i=d,s,(b) 

where we have defined 
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and s = s/rric, z; L = mi/m c . The logarithmic dependence on the internal quark mass mi in the 
second term of Equation 94 cancels against a similar term in the Inami-Lim function F\(xi) 
entering in Cg(Mjy), leaving no spurious divergences in the mi — > limit. 5 

5 Fajfer et al. (107) do not take the mixing of Og with O2 into account. This results in a prediction for the 
short-distance components that is mainly given by these logarithms. 
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3.1.1 THE c u£+£- DECAY RATES 



To compute the differential decay rate in terms of the Wilson coefficients, we use the two-loop 
QCD corrected value of C| ff (m c ) as obtained in Reference (98); we compute Cg ff (m c ) from 
Equation 94 and Cio(m c ) = C\q(Mw) from Equation 93. The differential decay rate in the 
approximation of massless leptons is given by 
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(96) 



where td refers to the lifetime of either D^ or D°. We estimate the inclusive branching ratios 
for m c = 1.5 GeV, m s = 0.15 GeV, m& = 4.8 GeV and = 0, 
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(97) 



The dominant contributions to the rates in Equation 97 come from the leading-order mixing of 
Og with the four-quark operators > the second term in Equation 94. When considering the 
contributions of various new-physics scenarios, one should remember that their magnitudes 
must be compared to the mixing of these operators. Shifts in the matching conditions for the 
Wilson coefficients CV, Cg, and Cio, even when large, may not be enough to give an observable 
deviation. 



3.1.2 THE c^u-f RATE 

The short-distance c — ► wy contribution to radiative charm decays was first studied in detail 
by Burdman et al. (97), who found that the effects of the leading logarithms on CV(mc) 
enhanced the branching ratio by several orders of magnitude. Even with such enhancement, 
the rates were too small. However, Greub et al. noted (98) that the leading logarithmic 
approximation was still affected by a fair amount of GIM suppression because the quark mass 
dependence on the resummed expressions was still mild. Going to two loops in the matrix 
elements of the operators in Equation 87, specifically in > leads to a more substantial mass 
dependence that in turn breaks GIM more efficiently. These two-loop contributions dominate 
the short-distance radiative amplitude, giving, for instance (98), 

B^ sd \D° -» X-y) ~ 2.5 x 10~ 8 . (98) 

Although this represents a very large enhancement even with respect to the leading logarithmic 
approximation (about five orders of magnitude!), it is still small, especially when compared 
with the estimated size of long-distance contributions (see below). 



3.1.3 EXCLUSIVE SEMILEPTONIC DECAY MODES 

The exclusive modes corresponding toc->u transitions are known to be dominated by long- 
distance dynamics. This is true for both the radiative and the semileptonic decays. For 
the D — > Xj exclusive modes, long-distance physics dominates all observables. However, in 
D — > Xl + l~ , it is in principle possible to escape the largest long-distance contributions by 
looking at regions of phase space away from resonances. We now discuss in some detail the 
computation of D — > tt£ + £~ and D — > p£ + £~ as presented in Reference (99). For completeness, 
the expectations in the exclusive radiative and neutrino modes are surveyed at the end of this 
section. 
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As a crude first estimate of the contributions of long-distance physics, we can consider the 
resonance process D — > XV — > X£ + £~, where V = 4>,p,uj. We isolate contributions from 
this particular mechanism by integrating dT/dq 2 over each resonance peak associated with 
an exchanged vector or pseudoscalar meson. The branching ratios thus obtained are in the 
0(1O~ 6 ) range (102). 

This result suggests that the long-distance contributions overwhelm the short-distance 
physics and any new physics that might be present. However, this is not always the case. 
A more thorough treatment requires looking at all the kinematically available regions in 
D — > X u £ + £~, not just the resonance region. The effect of these states can be thought of 
as a shift in the short-distance coefficient C| ff in Equation 94, since V — ► £ + £~ selects a 
vector coupling to the leptons. This follows from Reference (103), which incorporates in a 
similar manner the resonant contributions to b — > q£ + £~ decays via a dispersion relation for 
£ + £~ — > hadrons. This procedure is manifestly gauge- invariant. The new contribution can be 
written via the replacement (103) 

Of - Of + % £ Ki 2 mv ^ + *- , (99) 



or . m 



where the sum is over the various relevant resonances, my i and are the resonance mass 
and width, and the factor ki ~ 0(1) is a free parameter adjusted to fit the nonleptonic decays 
D — > XVi when the Vi are on shell. We obtain ~ 3.6, k p ~ 0.7, and ^ ~ 3.1. The 
latter result comes from assuming B(D + — > -k + uj) = 10 -3 , since a direct measurement is not 
available yet. 

D + — > n + e + e~: The main long-distance contributions come from the </>, p, and uj resonances. 
The 7] and rj' effects are negligible. The dilepton mass distribution for this decay takes the 
form (99) 



£=i£fti a i/ + Mi s 



+ |Ciof , ( 10 °) 



where s = m 2 e is the square of the dilepton mass. Here we have used the heavy-quark 
spin-symmetry relations that relate the matrix elements of Oj to the "semileptonic" matrix 
elements of Og and Oiq (104). An additional form factor is formally still present, but its 
contribution to the decay rate is suppressed by (m^/m/)) 2 and is neglected here. Precise 
measurements of D — ► -k£u will give us /+((/ 2 ). In the meantime, we make use of the prediction 
of chiral perturbation theory for heavy hadrons (ChPTHH) (105), which at low recoil gives 

/+( s = ~TT\ 7772~T- !01) 

f v (1 - s/M^) 

Here we use the recent CLEO measurement (106) gD*D-w = 0.59 ± 0.1 ± 0.07, and we take 
fo = 200 MeV. In Figure 7, we present this distribution as a function of the dilepton mass. 
The two narrow peaks are the <f> and the u>, which sits on top of the broader p. The total rate 
results in B(D + — ► 7r + e + e _ ) ~ 2 x 10~ 6 . Although most of this branching ratio arises from 
the intermediate tt + (/) state, we can see from Figure 7 that new-physics effects as low as 10 -7 
can be observed as long as such sensitivity is achieved in the regions away from the uj and 4> 
resonances, both at low and high dilepton mass squared. 

Fajfer et al. use a different approach to compute the long-distance effects in D — ► n£ + £~ 
(107). They estimate individual contributions to the amplitude by using a combination of 
vector meson dominance (VMD), factorization, ChPTHH, and hidden local symmetry. In 
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addition to the large inherent uncertainties associated with each of these methods, there is no 
clear prescription of how to fix the relative signs of the various contributions, nor — in contrast 
to the previous approach — the sign relative to the short-distance contribution. Within these 
large uncertainties, the results of Reference (107) and Reference (99) do not differ by much. 

D + — > p + e + e~: Once again, we follow closely the calculation of Burdman et al. (99). 
Because fewer data are currently available on the D — > VV' modes, we take the values of 
the Ki in Equation 99 from the fits to the D + — > ir + V case studied above. Again, once 
precise measurements of the D — > p£v form factors are available, heavy quark spin symmetry 
relations can be used to turn these into D — > p£ + £~ form factors. Lacking these at the 
moment, we use the extracted values from the D — > K*£u data (108) and assume SU(3) 
symmetry (113). The total integrated branching ratio is B(D° — > p°e + e~) = 1.8 x 10 -6 (i.e., 
B{D + — > p + e + e~) = 4.5 x 10~ 6 ). As Figure 8 shows, once again most of this rate comes 
from the resonance contributions. However, there is also a region — in this case confined to 
low values of m ee owing to the kinematics — where sensitive measurements could test the 
standard-model short-distance structure of these transitions. 

Fajfer et al. (109) studied these modes, following the approach of Reference (107). Their 
results (109) show a large enhancement at low m ee when compared with Figure 8. If this 
long-distance enhancement is present, it could dominate the low m ee region, rendering the 
D — > V£ + £~ modes useless to test short-distance physics at any level. Therefore, we examine 
this discrepancy carefully. A l/m 2 e enhancement can only appear as a result of nonfactorizable 
contributions. This is clear from References (110) and (111): the factorization amplitude for 
D — > pV, when combined with a gauge-invariant (7 — V) mixing, leads to a null contribution 
to D — > V£ + £~ . This is reflected also in the fact that the mixing of the operator O2 with 
Oj is nonfactorizable (111). A resonant contribution to O7, leading to a 1/q 2 behavior, 
is then proportional to Cf , which is mostly given by the O2 mixing. In addition, when 
compared with the usual short-distance matrix element of O7, this resonant contribution 
will be further suppressed by the factor gy (q 2 )A ni (q 2 ) , where gv{l 2 ) is the (7 — V) mixing 
form factor, and A ni (q 2 ) parameterizes the nonfactorizable amplitude (pV\Oj\D) , which is 
of C(Aqcd/"^c) (112). Thus, even if we take the on-shell values for these quantities, the 
resonant contribution to O7 is likely to be below 10% of the standard-model short-distance 
contribution. The actual off-shell values at low q 2 far from the resonances are likely to be even 
smaller. We then conclude that the 1/q 2 enhancement is mostly given by the short-distance 
contribution. This is only noticeable at extremely small values of the dilepton mass, so that 
it is likely to be beyond the experimental sensitivity in the electron modes (due to Dalitz 
conversion), whereas in the muon modes it lies beyond the physical region. On the other 
hand, the factorizable pieces contribute to the matrix elements of O9, just as in Equation 99, 
and give no enhancement at low values of q 2 . We then conclude that there should be no very 
large long-distance contributions at low m ee . 

The p modes also contain angular information in the form of a forward-backward asymmetry 
for the lepton pair. Because this asymmetry results from the interference between the vector 
and the axial- vector couplings of the leptons, it is negligible in the standard model, since 
vector couplings due to vector mesons overwhelm axial- vector couplings. This is true even 
away from the resonance region, since the coefficients Cg'^eff and C| ff get large enhancements 
due to mixing with O2 and QCD corrections, whereas C\q — the axial-vector coupling — is not 
affected by any of these, which results in a very small interference. 

For both the ir and p modes, the sensitivity to new-physics effects is reserved for large 0(1) 
enhancements because the long-distance contributions are still important even away from the 
resonances. Table 14 summarizes theoretical predictions from Reference (99). In addition, 
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some modes are almost exclusively long-distance physics. Examples are D° — » K°^£ + £ and 
the radiative D° — > K°^*'j, dominated by W exchange diagrams, as well as D ± -> K^t+f 
and the radiative — > 2f ± W'y, which contain both W annihilation and exchange. The 
measurements of these modes, although not directly constraining new physics, will help us 
understand long-distance physics. This may prove crucial to test the short-distance physics 
in the n and p modes. 



Exclusive Radiative Decays: Exclusive decays mediated by the c — ► wy transition are 
expected to be plagued by large hadronic uncertainties. As mentioned, the large mixing of 
the Oj operator with the four-quark operators, especially O2, and the propagation of light 
quarks in the loops indicate the presence of potentially large nonperturbative effects. These 
are not calculable from first principles nor in a controlled approximation (other than lattice 
gauge theory). Moreover, even if lattice computations of these effects become available, they 
typically overwhelm the standard-model short-distance contributions. Thus, modes such as 
D — > p7 are not expected to be a probe of the short-distance structure of the standard model 
to the extent B — > K*^f can be if the transition form factor is known precisely. 

On the other hand, one can try to estimate the size of the long-distance contributions and 
therefore the branching fractions of these modes. This is interesting in its own right; experi- 
mental observation of these modes will give us guidance in our otherwise limited understanding 
of these nonperturbative effects. 

Several attempts at estimating the long-distance contributions have been made (97, 115- 
118). An example is the decay D° — > p°7. We can identify two types of long-distance contri- 
butions: pole-mediated and vector-meson dominance (VMD) transitions. Pole contributions 
can be thought of as driven by "annihilation" diagrams. The effective weak Hamiltonian 
underlying these transitions is 

H w = Mul7X)(*'l7 m c l ) + a 2 (s' L ^d' L )(u L ^c L )] , (102) 

with d! = V u dd + V us s and s' = V cs s + V c dd. There are two pole diagrams giving (97) 

A(r>° _^ „°\ — ffp^T /„Qi it 1 rA 
A{L> — > p J pole — — g oK 11 \H W \D ) 

m B — 

+ 9 2 D *° D °l (p°\H w \D*°). (103) 
m D — m D * 

Here, H w refers to the four-fermion weak Hamiltonian governing nonleptonic weak decays. 
For instance, 

(AH W \D°) = -| fl2 K)»|H, (104) 

where a2(m c ) is determined in D nonleptonic decays. These two pole diagrams may can- 
cel rather efficiently. Since the details are extremely model-dependent, we could ask what 
branching ratio would be implied by one of them alone. This typically leads to (97) B(D° — > 
P°7)poie < few x 10~ 7 . However, QCD sum rules may be used (115) to compute the annihi- 
lation contributions, and these are found to give B(D° — ► p°j) ~ few x 10~ 6 . 

On the other hand, VMD contributions come from considering the nonleptonic intermediate 
states. In this case, this corresponds to D° — ► p°V — ► p 0r y, where the neutral vector boson 
V turns into an on-shell photon. Various methods have been used to compute the nonlep- 
tonic and V —> p amplitudes (97, 116, 117). A common assumption to estimate the VMD 
amplitude has been that of factorization (119). However, the contribution of the factorized 
nonleptonic amplitude vanishes when the photon is on-shell (110, 111). This is a consequence 
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of gauge invariance and is related to the fact that the mixing of four-quark operators with 
the photon penguin operator Oj vanishes unless nonfactorizable gluons are exchanged. Thus, 
nonfactorizable contributions to the nonleptonic amplitude constitute the leading effect in 
the VMD amplitude. It is therefore possible that the VMD contributions to weak radiative 
decays of charm mesons are overestimated. At the same time, it is possible that the charm 
quark is not heavy enough for the nonfactorizable effects to be suppressed. The suppression 
is formally 0(A/m c ), with A a typical scale of strong interactions. We conclude that uncer- 
tainties in these modes are still very large. Table 12 summarizes predictions for these decays 
in different model computations. The Belle collaboration recently announced a measurement 
of D° -> (fry, with a branching ratio of (138) B(D° 7 ) = (2.60±^;^±g;}f) x 10" 5 . As one 
can see from Table 11, this is consistent with the upper end of the predictions in Reference 
(97), which were obtained by making use of VMD plus the data on the relevant nonleptonic 
decay in addition to the pole contributions. If this trend is confirmed in this as well as other 
modes, this might point in the direction of large nonfactorizable contributions. Experimental 
bounds are closing in on some of these predictions and will undoubtedly shed light on the size 
of these long-distance effects. 

3.1.4 OTHER DECAY MODES 

Our discussion has neglected many decay modes to focus on those that can potentially test 
the short-distance structure of the standard model and are reachable in present or planned 
experiments. Here is a sample of other modes studied in the literature. 

jj° — > jry. The st andar d-model short-distance contributions can be obtained from the two- 
loop c — > amplitude. This results in (99) i? SD (D° — ► 77) ~ 3 x 10~ n . There are several 
types of long-range effects. Fajfer et al. (120) estimate these effects using ChPTHH to one 
loop. This gives B LD (D° -> 77) ~ (1 ± 0.5) x 10~ 8 . Burdman et al. (99) consider various 
long-distance effects and obtain similar results. In this case, the main contributions are found 
to come from VMD and the K + K~ unitarity contribution. 

D° — ► £ + £~: The short-distance contributions to this mode are also extremely suppressed, 
not only by helicity but also by the quark masses in the loop. Unlike in c — > wy, the mixing with 
O2 does not help. In Reference (99), the short distance is estimated at B (D°£ + £~) ^ 10~ 18 . 
The most important source of long-distance effects is the two-photon unitary contribution, 
which gives 

B hD (D° -> £ + r) ~ 3 x 10~ 5 B(D° -> 77). (105) 

D — > Xvv: Short- and long-distance contributions to c — > uvv processes in the standard 
model are extremely small, typically resulting in (99) B(c — ► uvv) 10~ 15 . 

3.2 Rare Charm Decays and New Physics 

As mentioned, charm-changing neutral-current processes such as D°-D° mixing and rare 
charm decays complement the constraints on extensions of the standard model obtained from 
processes initiated by down quarks, such as kaon and l?-meson transitions. Although we have 
seen that bounds on Amp are quite constraining in a variety of models, new physics may still 
show itself in rare charm decays. We mainly review the potential for signals from supersym- 
metric theories (with and without -R-parity conservation) and from new strong dynamics at 
the TeV scale. We briefly comment on the sensitivity to other new physics, such as theories 
with extra dimensions and extended gauge and matter sectors. 
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3.2.1 THE MINIMAL SUPERSYMMETRIC STANDARD MODEL 

The MSSM adds to the standard- model description of loop- mediated processes contributions 
due to gluino-squark exchange, chargino/neutralino-squark exchange, and charged Higgs- 
quark exchange. This last contribution carries the same CKM structure as in the standard- 
model loop diagram and is proportional to the internal and external quark masses; thus, its 
effects in rare charm transitions are small and we neglect it here. The gluino-squark con- 
tribution proceeds via flavor-diagonal vertices proportional to the strong coupling constant 
and in principle dominates the CKM-suppressed, weak-scale strength chargino/neutralino- 
squark contributions. We therefore consider only the case of gluino-squark exchange here as 
an estimate of the potential size of SUSY effects in rare charm decays. 

A typical squark-gluino contribution is depicted in Figure 9. The corresponding effects in 
the c — ► u transitions were studied for D — > X u ^ (17) and for D — ► X u £ + £~ (99). Within 
the context of the mass insertion approximation (14), the effects are included in the Wilson 
coefficients corresponding to the decay D — > X£ + l~ via 

d = Cf M + Cf , (106) 

for i = 7, 9, 10. Allowing for only one insertion, the explicit contributions from the gluino- 
squark diagrams are (17, 99, 121) 



n~a 8 \/2_ f Pi32(u) \ M §\ /in - 

° 7 = "9 G^Ml™ s |(^)"— T- + ^VrPMu) — ) , (107) 



and 

Cf = ~ ~^2^ (^ 2 )llP a2(u), (108) 

with the contribution to C\q vanishing at this order because of the helicity structure. 
If we allow for two mass insertions, there is a contribution to Cg^o given by 

C l = -\- (^)lr(S^)lrP 32(u) = (109) 
9 a 1 — 4 sin t)w 

Here, u = and the functions Pijk(u) are defined as 

pi T i(1 _ r )j 

P ijk (u) = / dx X[ X) (110) 

(1 — X + UX) K 







In addition, the operator basis can be extended by the "wrong chirality" operators O7, Og, and 
O10, obtained by switching the quark chiralities in Equations 90 and 91. The gluino-squark 
contributions to the corresponding Wilson coefficients are 

c? = -^^-^12)^^ (in) 

Cf = -^^2vra s (^ 2 )^042(n)-(l-4sin 2 MCf , 
In- 

Cf = -q- (^2)^2)1^032 

where the expression for C*f is again obtained with a double insertion. 

As noted in References (121) and (17), in both &} and Cy the term in which the squark 
chirality labels are mixed introduces the enhancement factor Mg/m c . In the standard model, 
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the chirality flip that appears in Oj occurs by a flip of one external quark line, resulting in 
a factor of m c included in the operator's definition. 6 However, in the gluino-squark diagram, 
the insertion of (Sf 2 )RL forces the chirality flip to take place in the gluino line, thus introduc- 
ing a Mg factor instead of m c . This yields a significant enhancement in the short-distance 
contributions to the process D — ► X u j (17), which is unfortunately obscured by the large 
long-range effects. This is not the case in c — > u£ + £~ processes. 

In order to estimate the effects in c — ► u£ + £~ transitions from the gluino contributions, we 
make use of the bounds of Table 1. Figures 10 and 11 show the dilepton mass distribution 
as a function of the dilepton mass for D + — ► 7r + e + e _ and D° — ► p°e + e~, respectively (from 
Reference (99)). Four sample cases are considered there (99): (I) M g = Mg = 250 GeV, (II) 
M g = 2Mg = 500 GeV, (III) M g = Mg = 1000 GeV, and (IV) M g = (1/2) Mg = 250 GeV. 

We first examine the D + — > ir + e + e~ case. Although the net effect is relatively small in 
the integrated rate (an increase ~ 20% or smaller), the enhancement due to the SUSY con- 
tributions is most conspicuous away from the vector resonances, particularly for low dilepton 
masses. Experiments sensitive to the dilepton mass distribution at the branching fraction of 
10 — 10 -8 can detect these SUSY contributions. However, the decays to a vector meson, 
such as D — > pe + e~, are more sensitive to the gluino exchange, as Figure 11 shows. 

The effect is quite pronounced and lies almost entirely in the low m ee region. This is mostly 
due to the contributions of {5\ 2 )rl to Cj and Cy, which contain the M g /m c enhancement as 
advertised above. This effect is intensified at low q 2 = m 2 e owing to the photon propagator 
(see, e.g., Equation 96 for the inclusive decays). This low-g 2 enhancement of the O-j contribu- 
tion is present in exclusive modes with vector mesons, such as D — > p£ + £~ , but not in modes 
with pseudoscalars, such as D — > n£ + £~ , since gauge invariance forces a cancellation of the 
l/q 2 factor in the latter case (see, e.g., Equation 100). This is apparent from a comparison 
of the low dilepton mass regions in Figures 10 and 11. Thus we see that rare charm decays 
are indeed sensitive to a generic extension of the standard model such as the MSSM. This is 
particularly true of the D — > p£ + £~ modes. 

3.2.2 SUPERS YMMETRY WITH i?-PARITY VIOLATION 

Imposing i?-parity conservation in the MSSM prohibits baryon- and lepton-number-violating 
terms in the superpotential. However, other symmetries can be invoked to prohibit rapid 
proton decay, such as baryon parity or lepton parity (see, e.g., (122)), and hence allow for 
-R-parity violation. The i?-parity-violating superpotential can be written as 



where L, Q, E, U, and D are the chiral superfields in the MSSM. The SU(3) color indices are 
denoted by a, (3, 7 = 1,2,3, the SU(2)^ indices by a, b = 1,2, and the generation indices by 
i,j, k = 1, 2, 3. The fields in Equation 112 are in the weak basis. 

The X' i j k term is the relevant one for the rare charm decays we consider here because it can 
give rise to tree-level squark-exchange contributions to decay channels such as D — ► X£ + £~, 
D — > £ + £~ , as well as the lepton-flavor-violating D — > Xp + e~ and D — ► p + e~ modes. Before 
considering the FCNC effects in D decays, we need to rotate the fields to the mass basis. This 
leads to 




(112) 



Wr p = X'ijk [NiVjiDt - EiUj] D k + -- 



(113) 



where V is the CKM matrix and we define 



A ijk — A irs ul rj ly sk • 



(114) 



The m u term, proportional to the (1 — 75) in the operator, is neglected. 
We ignore bilinear terms that are not relevant to our discussion of FCNC effects. 



T 
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Here, U L and V R are the matrices used to rotate the left-handed up- and right-handed down- 
quark fields to the mass basis. Written in terms of component fields, this interaction now 
reads 

wy = A^{^ 7 [444 + dl44 + (4r(4) c ^] 

-444 - 444 - (4)*(4) c 4} • (us) 

The last term in Equation 115 can give rise to the processes c — > u££^ at tree level via the 
exchange of a down squark. This leads to effects that are proportional to A^ 2fc A- lfc with i = 1, 2 
(owing to kinematical restrictions). 

Constraints on these coefficients have been derived (for a recent review, see (123)). For 
instance, Agashe &: Graesser (124) obtain tight bounds from K + — > -k^vv by assuming that 
only one i?-parity-violating coupling satisfies A^ fe 7^ 0. We update this bound by using 
the latest experimental result (125) B(K + -> tt+vv) = (l^^oi!) x 10~ 10 , which yields 
Kjk < 0-005. However, this bound can be avoided in the single-coupling scheme (124), where 
only one i?-parity-violating coupling is taken to be nonzero in the weak basis. In this case, 
it is possible that flavor rotations may restrict the i?-parity-breaking-induced flavor violation 
to be present in either the charge — | or +| quark sectors, but not both. Then large effects 
are possible in the up sector for observables such as D°-D° mixing and rare decays without 
affecting the down-quark sector. 

Agashe &: Graesser (124) obtain a rather loose constraint on the i?-parity-breaking cou- 
plings from D° mixing, which could result in large effects in c — > u££^ ) decays. Here, we take 
a conservative approach and use more model-independent bounds. The constraints on the 
i?-parity-breaking couplings for the processes of interest here are collected in Table 13 from 
Reference (123). The charged-current universality bounds assume three generations. The 
7r-decay constraint is given by the quantity R n = T{ir — > ev)/T(ir — > fiu). The limits obtained 
from D — > Kiv were first obtained in Reference (126). 

We first consider the contributions to c — > u£ + £~ . The tree- level exchange of down squarks 
results in the effective interaction (99) 

5H cS = -^kk (4)c C£ u L {£ L ) c , (116) 



which, after Fierz rearrangement, gives 



sn eS = (u L ^c L )(h^e L ). (ii7) 

m d R 

This corresponds to contributions at the high-energy scale to the Wilson coefficients Cg and 
C10 given by 

5C* = -5C 10 = S 4^(^) 2 K 2k K lk . (118) 



2a 2 



a R , 



If we now specify £ = e and use the bounds from Table 13, we arrive at the constraint 



5Ct = -5Cf < 1.10 (^] 



lk 



0.04 / \ 0.02 



(119) 



Notice that these are independent of the squark mass, which cancels. Taking this upper limit 
on the Wilson coefficients results in the dotted-dashed lines of Figures 7 and 8 corresponding 
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to D + — ► 7r + e + e _ and D° — ► p°e + e~ , respectively. The effect in these rates is small, of order 
10% at most, whereas the experimental bounds are a factor of 20 above this level in the most 
restrictive case (given by the pion mode). On the other hand, for £ = p,, we obtain 

iC -- iCfo < 17 . 4 (^)(||). (120) 

When the existing constraints on A' ljfc are taken into account, F-parity violation results in 

small deviations in D — > (ir, p)e + e~ . However, the bounds on X' 2 j k are loose and lead to very 
large effects in the £ = p modes. In fact, the allowed values from other observables saturate 
the current experimental limits for D — ► irp + p~ and D — ► pp + p~ , resulting in the bound (99) 

Am A' 21fc < 0.004. (121) 

These large effects would be clearly observable away from the resonances. 

In addition, the angular information in D —> pp + p~ can be used to confirm the new-physics 
origin of the large deviations. If we define the forward-backward asymmetry for leptons as 

A ^° dxdq 2 ^ X /-l dxdq 2 ^ X n 00 \ 

A-FB\q ) = 2t ' y 122 ) 

dq 2 

where x = cos#, with 9 being the angle between the £ + and the decaying D meson in the 
£ + £~ rest frame (as mentioned in Section 3.1.3), the standard-model prediction for this quan- 
tity is nearly zero. This is because the rate comes almost entirely from the long-distance 
contributions to the lepton vector coupling (e.g., Og), making the interference with the axial- 
vector couplings negligible. This can be seen by inspecting the numerator of the asymme- 
try (113, 127) 

A FB (q 2 ) ~4m D k do {Cf g f + ^Cf (f G - g F) j , (123) 

where h is the vector meson three- momentum in the D rest frame, and /, g, F, and G are 
form factors. Because the standard-model amplitude is dominated by the long-distance vector 
intermediate states, we have C| ff S> C\q. New-physics contributions that make C\q ~ Cg ff will 
hence generate a sizeable asymmetry. This is actually the case in i?-parity-violating SUSY. 
Burdman et al. (99) show that the asymmetry in D — > pp + p~ is predicted to be quite large 
for the allowed values of the -R-parity- violating parameters. 

The effective Wilson coefficients of Equation 118 also lead to a contribution to the two-body 
decay D° — > p + p~ . The imparity-violating contribution to the branching ratio then reads 



B^(D° - M V) = r D o flmlmn Jl- ^ ^ ^ ■ (124) 



2 



4m2 {^22kK 

Applying the bound in Equation 121 gives the constraint 8 

B **(D° n+n~) < 3.5 x 10- 7 • (125) 

Thus, i?-parity violation could give an effect that could soon be probed in these modes. 

8 In Reference (99), this expression (Equation 86) was incorrectly given. Also, the branching ratio stated 
there did not reflect the bound from Equation 122, but the less restrictive bounds to the individual couplings. 
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Next, we consider the products of imparity- violating couplings that lead to lepton flavor 
violation. For instance, the products A' llfc A2 2fc and X^ik^'uk w * n §^ ve r i se to D + — ► vr + /x + e _ . 
This leads to 



^f = -^ = 4.6x ^ ^ + ^ ^ , (126) 




which results in B^ V (D + — > 7r + /x + e _ ) < 3 x 10~ 5 . Here again, experiment is on the verge of 
being sensitive to i?-parity-violating effects. Similarly, for the corresponding two-body decay 
we have 



B ^(D° fj+e-) < 0.5 x 10~ 6 x M m + M M . (i 27 ) 




Last, we consider the i?-parity-violation contributions to — ► e^v and D^r — > e u. These 
decay modes receive large enhancements mostly from the s-channel exchange of sleptons (128). 
Unlike the t-channel squark exchange discussed above, this results in amplitudes that are 
unsuppressed by helicity. The largest lepton- number violation occurs in — > e^v T through 
the product | A231 A'221 1 , and in — > e^v T through | A231 A'222 1 • Here \jk is the strength of the 
cubic lepton superfield interactions in the i?-parity-violating superpotential term ea^L^L^E^. 
When current bounds are taken into account, it is found that i?-parity violation could result 
in (128) B^ p (D :iz — > e^v) ~ 10 -4 , several orders of magnitude above the standard-model 
prediction. Even more dramatically, we could have B^ p (D^ — > e^v) ~ 5 x 10~ 3 . Table 14 
summarizes the predictions from Reference (99) in both the standard model and i?-parity- 
violating SUSY. 

3.2.3 STRONG DYNAMICS 

Technicolor Models: In standard technicolor theories, both fermions and technifermions 
transform under the new gauge interaction of extended technicolor (ETC). As we saw in 
Section 2.3.2, this leads to the presence of four-quark operators coming from the diagonal 
ETC generators and characterized by a mass- scale M bounded by D° - D° mixing to be 
greater than 100 TeV or so. However, additional operators are generated at low energies 
that are not suppressed by M. The condensation of technifermions leading to electroweak 
symmetry breaking leads to fermion mass terms of the form 

m q c^{fT) TC . (128) 

iw ETC 

Operators arising from the technifermion interactions have been shown (129) to give rise to 
FCNC involving the Z boson, 

^-U?Z» {u Ll ,c L ) and C 2 ^— ^-U^Z^ (u Ll ,c L ), (129) 
87rv sm 2vw 8nv sin 2f)y/ 

where Ul is the unitary matrix rotating left-handed up-type quark fields into their mass basis 
and £ is a model-dependent quantity of 0(1). The induced flavor-conserving Z coupling was 
first studied in Reference (129), and flavor-changing effects in B decays have been examined 
in References (130) and (131). 

The flavor-changing vertices in Equation 129 induce contributions to c — > u£ + £~. These 
appear mostly as a shift in the Wilson coefficient Cio(Mw), 

2v a 
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where we assume Uf^ 1 ~ A ~ 0.22 (i.e., one power of the Cabibbo angle) and we take m c = 
1.4 GeV. Although this represents a very large enhancement with respect to the standard- 
model value of C\o(M\y), it does not translate into a large deviation in the branching ratio. 
As we have seen, these are dominated by the mixing of the operator O2 with O9, leading to 
a very large value of Cg ff . The contribution in Equation 130 represents only a few-percent 
effect in the branching ratio with respect to the standard model. However, the effect is quite 
large in the region of low dilepton mass. 

Furthermore, the interaction in Equation 129 can also mediate D° — > fj, + fi". The corre- 
sponding amplitude is 



•4dV>- - U cu |^ -j= sin 2 9 w f D . (131) 



This results in the branching ratio I? ETC (D — ► /x + /i~) ~ 0.6 x 10~ 10 , which, although 
still small, is not only several orders of magnitude larger than the standard-model short- 
distance contribution but also more than two orders of magnitude larger than the long-distance 
estimates. 

Finally, the FCNC vertices of the Z boson in Equation 129 also give large contributions to 
c — ► uvv. The enhancement is considerable and results in the branching ratio 

^ C ^p-e 4 (^f) 2 2xl0^. (132) 

Topcolor: As discussed in Section 2.3.2, the topcolor interactions must be nonuniversal, 
and so they mediate FCNCs at tree level. After the rotations of quark fields defined in 
Equation 40, the exchange of top-gluons generates four-fermion couplings of the form 

47rQ ^ C 2 ° tg2 U tc *U tu (u lfl T a t)(t^T a c) (133) 

G 

where U lJ = U 1 ^ + IA 1 ^ and M is the mass of the exchanged color-octet gauge boson. The 
one-loop insertion of this operator results in contributions to the operators Og and Oiq in 
c — > u£ + £~ as well as in the purely leptonic decays. These could lead to large effects when 
compared to the standard model 



sin 5 9 C 



/lTeVV , s 

Us-) ' (134) 



where the effects are rather modest unless the quark rotation matrices are larger than ex- 
pected. This would not be unnatural, for instance, for Ur, since the rotations of right-handed 
quarks are not related to any known observable in the standard model. 



3.2.4 OTHER NEW-PHYSICS SCENARIOS 

Extensions of the standard model leading to effects in rare charm decays also tend to result 
in large contributions to D° — D° mixing. Burdman et al. have evaluated a long list of 
these scenarios in detail (99). Generically, the effects are either negligible or amount to O(l) 
enhancement over the standard-model short-distance contributions. 

As mentioned in Section 2.3.3, compact extra dimensions may lead not only to massive 
scalar and fermionic states but also to nonuniversal couplings of the standard-model fermions 
to Kaluza-Klein excitations of gauge bosons that may induce flavor-violating loop effects. In 
general, the largest effects in rare charm decays are associated with massive neutral gauge 
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bosons (e.g., KK excitation of a Z'), which generate a FCNC current in the up-quark sector 
and then decay into either charged leptons or neutrinos. With masses starting around the TeV 
scale, these states could lead to O(l) enhancements in c — ► u£ + £~ modes, when compared to 
the standard-model short-distance predictions. Thus, in the charged-lepton modes this can 
translate into observable effects in the low rri£g window. The enhancement in the c — > uvv 
modes could be several orders of magnitude above the standard-model predictions, although 
they may be very difficult to observe. 

Many other new-physics scenarios involving additional matter and/or gauge fields would 
lead to contributions to flavor physics and in particular to rare charm decays. Most contribu- 
tions that are potentially large tend to be constrained by D° — D° mixing. Such is the case, 
for instance, with extra down-type quarks and extra gauge bosons. 

3.3 Experimental Bounds and Prospects 

Experiments that can measure rare D decay branching ratios at the level of a part per million 
will start to confront models of new physics in an interesting way. Although current experi- 
mental limits are, in most cases, well above this level of sensitivity, the arsenal of experiments 
opening up from CDF to the B factories and charm factory may allow observation of rare 
decays within the decade. 

We first review the present experimental status of charm rare decays. We then provide 
estimates of the reach of the currently running experiments CLEO-c and BaBar /Belle made 
by the CLEO-c collaboration, estimate the rare-decay reach of future facilities, and finally 
briefly consider the experimental dilepton mass resolution in decays of the type D — > ire + e~ . 

Fixed-target experiments have traditionally had the best rare-decay limits for low-multiplicity 
final states involving charged particles. In fixed-target experiments, rare decays are isolated 
using the same procedure used to isolate nonrare charm decays. The procedure is described 
in Section 2.4.6 of this review. 

Table 15 and Figure 12 show some of the best recent results and compare them to theory. 
At CDF, the search strategy for rare charm decays is still being refined but it will be similar to 
that developed for rare B decays such as B + — ► K + ^ + fi~ (87). Recently CDF has produced 
the best upper limit for D° — > u + u~ (88). 

At e + e~ B factories operating at the T(45), the background to rare charm decays is very 
large, and so far searches have been performed only for D° modes where the D* tag can 
be used to remove the significant combinatoric background. Using this technique, the first 
observation of a FCNC mode in the charm system has recently been made by BELLE (138), as 
shown in Figure 13. The measured branching ratio is B(D° — ► <jyy) = (2.60^ogi'toi7) x ^ 

Tables 16 and 17 list a collection of rare D° nodes. The second column is the best limit 
currently available from an e + e~ experiment running at or near the T(4S) resonance. The 
third column shows cases where the best limit was not obtained on the T(45). The next-to-last 
column estimates the sensitivity of an e + e~ experiment with 400 fb _1 of integrated luminosity 
at or near the T(45), based on the efficiencies and background levels found by CLEO, and the 
final column shows the projected CLEO-c limits based on Monte Carlo estimates of efficiency 
and tagging rates, and assuming an integrated luminosity of 3 fb _1 on the ^(3770) with a 10% 
tagging probability. This amount of data could be accumulated in one year. At the ^(3770) 
the assumption is made for all modes that the signal is background-free. 

Tables 18 and 19 show the current experimental status as well as projected sensitivities 
for a collection of rare D + decays. In these cases, there are no recent limits available from 
e + e~ experiments, so only the current world best measurement and the projected CLEO-c 
sensitivity [again based on 3 fb _1 of integrated luminosity on the ^(3770)] are shown. 

If tagging is not needed, the limits at CLEO-c will improve by a factor of ten. Although 
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these estimated limits are quite rough, we can draw some general conclusions from them. In 
general, CLEO-c and the £>-factory experiments have complementary strengths. For some 
final states, CDF will have the greatest sensitivity. 9 

The most reliable way to compare the reach of future experiments for rare decays would be 
to compile a list of sensitivities estimated by each experiment based on detailed simulations. 
However, as with 5Acp, in almost all cases the simulations or other detailed estimates are 
nonexistent. Instead, to compare the reach of future experiments for rare decays we again 
use the number of D° — > K~ir + reconstructed by each experiment. The mode we choose for 
the extrapolation for hadron experiments is D + — > n + p + p~ . The normalization is provided 
by the E791 search for that mode, which found B{D + -> 7r + //+/x~) < 5.2 x 1C T 5 at 90% CL. 
Because background is present, we scale this result by 1/ \/N(D^ — I K~tt + ). For the e + e~ 
experiments, the mode we choose for the extrapolation is D° — > n°e + e~. The normalization 
for machines operating in the U-meson threshold region is the result of the CLEO II search, 
which found B(D° -> vr°e + e-) < 4.5 x 10" 5 at 90% CL. We scale this by l/ v /iV(L' u K~ir+) 
in order to obtain upper limits for the £>-factory and super-i?-factory experiments. For the 
charm factories, reconstructing the second charm meson in the event ensures a background-free 
situation and the scaling is 1 /L using the CLEO-c estimated sensitivity as the normalization. 
The results of this exercise are in Table 11. 

To estimate the number of events each experiment will observe, we use B(D + — ► 7r + /i + /i _ ) ~ 
6.2 x 10 -6 and B(D° — > 7r°e + e~) ~ 8 x 10~ 7 . As can be seen from the results in Table 11, 
sensitivity to interesting physics in these modes can be achieved by BES III, BTeV, or the 
"super" machines. 

For the case of the dilepton modes D — > tt£ + £~ and D — > p£ + £~, as emphasized in Sec- 
tion 3.1.3, the sensitivity to short-distance physics, whether it is that of the standard model or 
one of its extensions, is in the region of low dilepton mass away from the resonances. Figure 14 
shows our Monte Carlo simulation of the experimental dilepton mass resolution in the decay 
D + — > ir + e + e~ for D mesons produced at the T(45) and for D mesons produced at thresh- 
old. The dilepton mass resolution is excellent in both cases. The dilepton mass resolution 
will also be very good at hadron facilities. In all cases, the dilepton resolution is expected to 
be adequate to sort between the standard-model and new-physics predictions for this class of 
decays. 

4 SUMMARY 

D°~D° mixing and rare charm decays will be further probed by experiment in the coming 
years. We have examined the potential of these measurements to test the standard model and 
its extensions. In the case of D°-D° mixing, the irreducible background from long-distance 
physics makes it difficult for a future positive observation to be attributed to short-distance 
physics. However, tight limits on Amp result in very important complementary constraints 
on physics beyond the standard model. On the other hand, although they are also generally 
affected by large theoretical uncertainties, we have seen that some rare charm decay modes 
can be used to probe short-distance physics. Such is the case with c — ► u£ + £~ decay modes. 
Sensitivity to the nonresonant regions in these modes could be achieved by experiments in 
the near future. In sum, these charm physics observables, when considered together with the 
results from the various B and kaon physics experiments, constitute a necessary low-energy 

9 We note that there is no commonly agreed prescription for calculating expected experimental upper limits. 
For an alternative, and more optimistic, estimate of the sensitivity to rare charm decays at a B factory, see 
Reference (57). 
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complement to direct searches for new physics at high energies. 
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Figure 1: Long-distance contribution from two charged pseudoscalar intermediate states. 
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Figure 2: Predictions for D° — D° mixing in the standard model and in some new-physics 
scenarios. From Reference (13). A: Standard-model predictions for x. □: standard-model 
predictions for y. •: New-physics predictions for x. The horizontal axis corresponds to the 
references in Reference (13). 
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Table 1: Bounds on (S^ll, {5i 2 )lr from D°-D° mixing (17) (neglecting the strong phase). 
All constraints should be multiplied by (Mg/500 GeV). 



M|/M| 






0.3 


0.03 


0.04 


1.0 


0.06 


0.02 


4.0 


0.14 


0.02 



Table 2: A comparison of the LEP experiments (35), CDF (36), E791 (37), FOCUS 
(38),CLEO (39), BaBar (40), and Belle (41). Each experiment has performed searches 
for charm mixing, rare decays, and CP violation. K~tt + is the number of reconstructed 
D° — > K~ir + used in published measurements. The proper time resolution for charm hadrons 
is denoted in at- 





Fixed Target 


e + e 


PP 




E791 


FOCUS 


LEP 


CLEO 


BaBar/Belle 


CDF 


Beam 


Hadron 


Photon 


e+e" - Z u 


e 


+e" 


PP 


K-7T+ 


~ 2 x 10 4 


~ 2 x 10 s 


~ 10 4 /expt. 


~ 2 x 10 & 


~ 10 b 


~ 5 x 10 b 




~ 40 fs 


~40 fs 


~ 100 fs 


~ 140 fs 


~ 160 fs 


~ 50 fs 



Table 3: Expected charm data sets for existing experiments and proposed facilities. "Current" 
is the size of the accumulated data set used in published physics analyses. "Full" is the size 
of the proposed complete data test, or the quantity of data collected per 10 s as indicated. 
Super Charm is not a proposed facility but is included for comparison. K + ir~ is the number 
of reconstructed D° — > K~ir + in the full data set. 



Experiment 


Current 


Full 


K--K+ 


BaBar 


91 fb" 1 


500 fb" 1 


6.6 x 10 b 


Belle 


46.2 fbr 1 


500 ftr 1 


6.6 x 10 b 


CDF (Run II) 


65 pb _i 


4.4 fb" 1 


30 x 10 b 


CLEO-c 




3 fb" 1 


5.5 x 10 b 


BESIII 




30 fb" 1 


5.5 x 10 b 


Super Charm 




500 ftr 1 


9.2 x 10710' s 


SuperKEKB 




2 ab" 1 


2.5 x 10 Y /10 Y s 


SuperBaBar 




10 ab^ 1 


1.3 x 10710' s 


BTeV 






~ 6 x 10 8 /10' s 
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Table 4: Comparison of the 95%-CL limits for the fit output parameters of BaBar, CLEO, and FOCUS 
when CP conservation is assumed in the fit, and when CP violation is allowed (denoted by CPV). The 
FOCUS entries are one dimensional limits. BaBar's limits include systematic uncertainties and were 
obtained in a fit that allowed x' 2 < 0. 





Rd 


y' 


x' 


x ri 


CLEO [%] 


(0.24, 0.69) 


((-5.2,0.2) 


(-2.8,2.8) 


< 0.0076 


BaBar [%] 


(0.24, 0.49) 


((-2.7,2.2) 




< 0.20 


FOCUS [%] 


(0.43, 1.78) 


(-12.4, -0.6) 


(-3.9,3.9) 




CLEO (CPV) [%] 


(0.24, 0.71) 


' (-5.8,1) ' 


(-2.9,2.9) 


< 0.082 


BaBar (CPV) [%] 


(0.23, 0.52) 


(-5.6,3.9) 




< 0.22 



Table 5: Measurements of the time-integrated rate for the wrong-sign decay D — ► K + ir , 
normalized to the right sign decay D° — ► K~tt + 





K-TT+ 


K + 7T~ 


Rws [%] 


A D [%] 


E791 (66) 


5.6K 


not quoted 


0.68+^ ± 0.07 




ALEPH (67) 


1038 


19 


1.84 ±0.59 ±0.34 




FOCUS (68) 


37K 


150 


0.404 ±0.085 ±0.025 




CLEO (61) 


13.5K 


45 


0.332t£o65 ± 0.040 


-2+ 19 ± 1 


Belle (63) 


83K 


317 


0.372 ± 0.025t£oi4 




BaBar (62) 


120K 


430 


0.357 ±0.022 ±0.027 


9.5 ±6.1 ±8.3 


Average 






0.368 ±0.021 





Table 6: Current D° lifetime difference measurements. BaBar (New) is a measurement of 
y cos <j) and supersedes BaBar. Only the former is used to compute the world average value of 
VCP- 





K~n + 


K-K+ 


7T 7T + 


VCP 


ref. 


E791 


35K 


3.2K 




(0.8 ± 2.9 ± 1.0)% 


(76) 


FOCUS 


120K 


10K 




(3.4 ± 1.4 ±0.7)% 


(77) 


CLEO 


20K 


1.9K 


0.7K 


(-1.1 ±2.5 ± 1.4)% 


(78) 


BaBar 


158K 


16. 5K 


8.4K 


(1.4 ±1.0+^)% 


(80) 


Belle 


214K 


18.3K 




(-0.5 ± 1.0 ±0.8)% 


(79) 


BaBar (New) 


220K 


26K 


12. 8K 


(0.8 ± 0.4±+^)% 


(81) 



Table 7: Comparison of measurements in Acp for D° modes, from E791 (89), FO- 
CUS (90), CDF (88), and CLEO (78) 





Mode 


A C p 


Mode 


A C p 


CLEO 


L> u -» K+K- 


(0.0 ±2.2 ±0.8)% 




(1.9 ±3.2 ±0.8)% 


E791 


L> u -> K+K- 


(-1.0 ± 4.9 ± 1.2)% 


D u -» vr+TT- 


(-4.9 ±7.8 ±2.5)% 


FOCUS 


D u K+K- 


(-0.1 ±2.2 ±1.5)% 


D U 7T+TT- 


(4.8 ±3.9 ±2.5)% 


CDF 


D u K+K~ 


(2.0 ± 1.7 ±0.6)% 




(3.0 ± 1.9 ±0.6)% 


CLEO 




(0.1 ±1.3)% 




(0.1 ±4.8)% 


CLEO 


D" - K S K» S 


(-23 ± 19)% 
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Figure 3: Our present knowledge of D°-D° mixing. The solid vertical lines indicate a "typical" 
standard-model prediction for x. The dashed vertical lines indicate the upper range of non- 
standard-model predictions for x. The horizontal band is the world average 95% CL limit in 
y. The circle with horizontal shading is the 95% CL limit in (x,y). The strong-phase shift 
Sktt between the Cabibbo-favored and doubly Cabibbo-suppressed decays is assumed to be 
zero in plotting the D° — > K + tt~ results, where in each case CP conservation is assumed. 
For the CLEO and FOCUS measurements, the statistical error is included; for the BaBar 
measurement, both the statistical and systematic errors are included. The strong phase shift 
is expected to be close to zero, but until it is actually measured, the allowed region from the 
D° — > K + ir~ measurements must be expanded to include the area swept out by rotating these 
regions about the origin. The three circles (small radius dashed, dot-dashed, and large radius 
dashed) are 2ir rotations of the BABAR, CLEO, and FOCUS regions, respectively. 
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Figure 4: Comparison of experimental determinations of yep- The average value is indicated 
by the shaded band. We calculate the average of the measurements to be (yep) = (0.8±0.5)%. 



Table 8: Comparison of measurements in Aqp for D + modes, from E791 (89) and FOCUS (90) 





Mode 


A C p 


FOCUS 


D+ K~K + ir + 


(0.6 ±1.1 ±0.5)% 




D+ -> K s tt + 


(-1.6 ± 1.5 ±0.9)% 




D+ -» K S K + w.r.t. (A^7r + 7r+) 


(6.9 ± 6.0 ± 1.5)% 




D+ K S K+ w.r.t. (K S 7T + ) 


(7.1 ±6.1 ±1.2)% 


E791 


D+ K-K+7T+ 


(-1.4 ±2.9)% 




D+ 4>ir + 


(-2.8 ±3.6)% 




D+ _^ K* U (892)K + 


(-1.0 ±5.0)% 




D + — > ir + ir~ir + 


(-1.7 ±4.2)% 



Table 9: Yields for direct CP violation in ^(3770) -> (CP Eigenstate l)(CP Eigenstate 2) 
in one year at CLEO-c (52) 



CP Eigenstate 1 


CP Eigenstate 2 


Events for 100% CP 


K+K- 


K + K- 


174 


K+K- 


A° 


171 


K s ir° 




183 


K S K S 


K+K- 


136 
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Figure 5: Measurements of Aqp for decay modes of the D . 



Table 10: Direct CP-violation asymmetry 5A reach from ^(3770) — > (flavor tag) {CP eigen- 
state), in one year at CLEO-c (52) 



CP Eigenstate 


Flavor- Tagged Sample 


5A 


K + K~ 


10200 


0.01 


K S TT° 


10400 


0.01 


K s u 


3500 


0.02 



54 



BURDMAN & SHIPSEY 





15 






a. 


10 


u 




< 


5 









-5 




-10 




-15 



6 7 8 

Mode Number 



1. FOCUS (D + -» K + K" x+) 

2. FOCUS (D + -» K s 71) 

3. FOCUS (D + -> K s K + (K TtV)) 

4. FOCUS (D + -> K s K + (K s n )) 

5. E791 (D + ^K K 7t ) 

6. E791 (D + ^(|)7t + ) 

7. E791 (D + -> K*(892) K + ) 

8. E791 (D + -> ji + it + 7I~ ) 



Figure 6: Measurements of Aqp for decay modes of the D 



Table 11: Projected CP asymmetry and sensitivity to rare decays for the full data set of 
existing experiments. For next-generation experiments, the sensitivities and rates are per year 
of 10 7 s. The yields at charm factories are computed without requiring a charm tag. Super 
Charm is not a proposed facility but is included for comparison. All entries are estimates by 
the authors. 



Experiment 


SA CP 


D -► tt£ + £~ limit 


D ire + £- yield 


BaBar 


3 x 1(T 3 


1 x 10~ b 


6 


Belle 


3 x 1(T 3 


1 x 10~ b 


6 


CDF (Run II) 


1 x l(T a 


2 x 10" b 


4 


CLEO-c 


0.01 


4 x 10~ b 


4 


BES III 


3 x 10" 3 


4 x 10~ Y 


40 


Super Charm 


2 x 10~ 4 


2 x 10" s 


700 


Super-KEK-B 


1 x W- ]i 


5 x 10-' 


60 


Super BaBar 


6 x 10~ 4 


3 x 10~ Y 


300 


BTeV 


3 x 10- 4 


4 x 10~ Y 


200 
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D*->tt* e* e 




m„ [GeV] 

Figure 7: The dilepton mass distribution for the D + — > n + e + e~, normalized to T D +. The 
solid line shows the sum of the short- and long-distance standard-model contributions. The 
dashed line represents the short-distance contribution only. The dotted-dashed line includes 
the contribution of ii-parity- violating terms in SUSY (see Section 3.2.2). 




m„ [GeV] 

Figure 8: The dilepton mass distribution for the D° — ► p°e + e~, normalized to T^o. The 
solid line shows the sum of the short- and long-distance standard-model contributions. The 
dashed line represents the short-distance contribution only. The dotted-dashed line includes 
the contribution of ii-parity- violating terms in SUSY (see Section 3.2.2). 



Table 12: Theoretical predictions for charm radiative decays (Br x 10 6 ) 



D —* Vj 


Reference (97) 


Reference (117) 


Reference (115) 


Df - p+ 1 


6-38 


20-80 


4.4 


D° -> K*°j 


7-12 


6-36 


0.18 


D° p° 7 


0.1 - 0.5 


0.1 - 1 


0.38 


D° w°7 


~ 0.2 


0.1 - 0.9 




D° (jp-i 


0.1 - 3.4 


0.4-0.9 




D+ -» p+7 


2-6 


0.4-6.3 


0.43 


D+ K* + j 


0.8 - 3 


1.2 - 5.1 
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Figure 9: A typical contribution to c — > u FCNC transitions in the MSSM. The cross denotes 
one mass insertion (5 12 )\\> , with A, A' = L, R. 



D*->tt* e* e 




m„ [GeV] 



Figure 10: The dilepton mass distribution for D + — > ir + e + e (normalized to T D +) in the 
MSSM with nonuniversal soft breaking effects. The solid line is the standard model. (I) 
M~ g = Mg = 250 GeV; (II) M~ g = 1M~ q = 500 GeV; (III) M g = M- q = 1000 GeV; (IV) 
M g = (1/2) Mq = 250 GeV. 



Table 13: Most stringent (2a) bounds for the i?-parity violation couplings entering in rare D 
decays, from (a) charged-current universality, (b) R n , and (c) D — > Klv* . See Ref. (123) for 
details. All numbers should be multiplied by (mi k /TOO GeV). 









A 22k 


0.02^ 


0.04( a ) 


0.06( b ) 


0.21^ 
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Figure 11: The dilepton mass distribution for D° — ► p°e + e (normalized to T D o) in the 
MSSM with nonuniversal soft breaking effects. The solid line is the standard model. (I) 
Mg = Mg = 250 GeV; (II) M g = 2M- q = 500 GeV; (III) M g = Mg = 1000 GeV; (IV) 
M g = (1/2) Mg = 250 GeV. 



Table 14: Comparison of various decay modes between the standard model and i?-parity 
violation. The third column shows how large the ii-parity-violating effect can be. 



Decay Mode 


Standard Model 


ftp 


D+ -> vr+e+e" 


2.0 x 10" 13 


2.3 x 10~ e 


D u -> p u e + e" 


1.8 x 10-° 


5.1 x 10~ b 


D+ -► n + fi + fl~ 


1.9 x lO" 13 


1.5 x 10" b 


D [ > -> pV + M" 


1.8 x 10 _B 


8.7 x 10~ e 


D u -> 


3.0 x 10~ i3 


3.5 x 10~ Y 


L> u -> e+e~ 




1.0 x 10~ iU 


D u -> n+e' 





1.0 x 10~ b 


D+ -> 7r+/i+e- 





3.0 x 10" & 


L> u -► pV+e" 





1.4 x 10" & 



Table 15: FOCUS results with incorporated systematic errors for the modes shown (from 
Reference (140)) Each number represents a 90% upper limit for the branching ratio of the 
decay mode listed. 



Decay Mode 


Result 


D+- 




9.2 x 10' 


-b 


D+ - 


■> K~Lt + [i + 


13.1 x 10 


-6 


D+ - 




8.8 x 10' 


■6 


D+ - 




4.8 x 10' 


6 


D + s - 


* K+fi+fj,- 


3.6 x 10- 


■5 


D + 3- 




1.3 x 10" 


-5 


D + s~ 




2.6 x 10- 


■5 




-4 7T~/i + /i + 


2.9 x 10" 


-5 
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Figure 12: FOCUS limits and sensitivity from Reference (140), compared to other experi- 
ments and theory. (For a treatment of sensitivity, see (139).) The previous limits, except 
for the E687 D + — > K~fi + fi + (141), are from the Fermilab experiment E791 (133). The 
theory estimates come from Reference (99) (i?-parity MSSM and standard modeUl) and Ref- 
erence (107) (standard model-2). The results for four-body decays measured by E791 (142) 
are not shown. 



Table 16: Current and projected 90%-CL upper limits in parts per million on rare D° decay 
modes. 



Mode 


T(4S)CLEO 


Best Upper Limit 


B Factory 


CLEO-c 




Reference 




Reference 








e + e~ 


(132) 


13 


(133) 


6.2 


0.6 


1.2 


H + fi~ 


(132) 


34 


(88) 


2.4 


3.4 


1.2 




(132) 


19 


(133) 


8.1 


1.9 


1.2 


■K°e + e~ 


(132) 


45 






1.0 


4 


vrV/i- 


(132) 


540 


(135) 


180 


53 


4 




(132) 


86 






8.5 


4 


K°e+e~ 


(132) 


110 






2.6 


7 




(132) 


670 


(135) 


260 


46 


7 




(132) 


100 






2.3 


7 


rje + e~ 


(132) 


110 






1.1 


10 


r]fi + fi- 


(132) 


530 






5.1 


10 




(132) 


100 






1 


10 
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Table 17: Current and projected 90%-CL upper limits in parts per million on rare D 1 
modes, (with the exception of (jyy, where observation has been recently reported) 





IV' . ' I ' 1 1' 1 




KOCT 1 T A T ^ QT" 1,1 V ) 1 n T" 




PT.F.O r 

V 1 j 1 j V J \ 




pe e 


( 1 10\ 

(13/ J 


1UU 




Q Q 


Z.O 


PI 1 I 1 








48 
*±o 




pfi e T 


{J-OZ) 








9 ^ 


OJG 6 


( 1 ^9^ 


lOU 




1 7 


Q 




(132) 


830 




57 


9 


f i) 1 1 p~F 

UJ Lb c 


(132) 


120 




11 


q 


K*(892)°e + e" 


(132) 


140 


(136) 24 


14 


4 


K*(892)°p + p- 


(132) 


1180 


(136) 47 


54 


4 


2T s (892)V ± e T 


(132) 


100 


(136) 83 


6.9 


4 


<pe + e~ 


(132) 


52 




5.1 


5 


(j)p + pT 


(132) 


410 


(136) 31 


28 


5 




(132) 


34 




3.4 


5 


7r + 7r _ 7r° p + fi~ 


none 




(135) 810 




6 


PI 


(137) 


240 




26 


1.5 




(137) 


240 




26 


8 




(138) 26i£l; 5 7 




21 


4 


A'*(892)°7 


(137) 


760 




380 


12 


77 


(134) 


29 






1.2 



Table 18: Current and projected 90%-CL upper limits on rare D + decay modes. 



Mode 


Best Upper Limit 
Reference 


CLEO-c 


7r + e + e~ 


(133) 


52 


1.5 




(140) 


8.8 


1.5 




(133) 


34 


1.5 


ir~e + e~ 


(133) 


96 


1.5 


ir~p + p + 


(140) 


4.8 


1.5 


tt~ p + e + 


(133) 


50 


1.5 


K+e+e- 


(133) 


200 


1.5 


K+fi+p- 


(140) 


9.2 


1.5 


K + p ± e T 


(133) 


68 


1.5 


K~e+e + 


(141) 


120 


1.5 


K~p + p + 


(140) 


13.1 


1.5 


K~p + e + 


(141) 


130 


1.5 
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Table 19: Current and projected 90%-CL upper limits in parts per million on rare D + decay 
modes. 



Mode 


Best Upper Limit 
Reference 


CLEO-c 


p + e + e~ 




7 




(135) 560 


7 


p + p^e^ 




7 


p~e + e + 




7 


p~p + p + 


(135) 560 


7 


p~ p + e + 




7 


K*(892)+e + e- 




10 


K*(892)+p+p- 




10 


K*(892)+p ± e T 




10 


K*(892)-e+e+ 




10 


K* (892)" /i+n+ 


(135) 850 


10 


K*{892)-p+e + 




10 
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Figure 13: The observation of D° ->■ 07, <f> -> K + K~ by the BELLE collaboration ((138)). 
(a) The invariant mass distribution of </>7 pairs for data (points with error bars), the fit to the 
data (line), the background component of the fit (dashed line), (fnr° background (dark shaded 
histogram), and the sum of (jyq and D + — ► 7r + 7r° backgrounds (light shaded histogram), (b) 
The helicity angle Ohel distribution in the signal region, where Ohel is defined as the angle 
between the direction of the K + momentum vector and the direction of the D° momentum 
vector in the rest frame of the (j) meson. Data (points with error bars), MC prediction (line), 
total background (light shaded histogram), and non (j>n background (dark shaded histogram, 
(c) Background subtracted 9h e i distribution in the signal region (points with error bars). The 
MC prediction is the line. The data are consistent with a sin 2 O^ei distribution as expected 
for D° — > (fry. 
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Figure 14: Dilepton mass resolution in D + — > ir + e + e for Ds produced at the T(4S) (solid) 
and at the "0(3770) (dashed). 



